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Even though matrix model partition functions do not exhaust the entire set of r-functions 
relevant for string theory, they seem to be elementary building blocks for many others and they 
seem to properly capture the fundamental symplicial nature of quantum gravity and string theory. 
We propose to consider matrix model partition functions as new special functions. This means they 
should be investigated and put into some standard form, with no reference to particular applications. 
At the same time, the tables and lists of properties should be full enough to avoid discoveries of 
unexpected peculiarities in new applications. This is a big job, and the present paper is just a step 
in this direction. Here we restrict our consideration to the finite-size Hermitean 1-matrix model and 
concentrate mostly on its phase/branch structure arising when the partition function is considered 
as a D-module. We discuss the role of the CIV-DV prepotential (as generating a possible basis in 
the linear space of solutions to the Virasoro constraints, but with a lack of understanding of why 
and how this basis is distinguished) and evaluate first few multiloop correlators, which generalize 
semicircular distribution to the case of multitrace and non-planar correlators. 
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Part I 

A Piece of Theory: Partition Functions and 
their Properties 

1 Introduction 

• One of the goals of generic string theory \1[ is to identify the properties of partition functions 
of various string models. These are defined as generating functions of all the correlators in a given 
quantum theory. Associated with the three existing formulations of quantum mechanics (linear algebra 
of operators in the Hilbert and Fock spaces, wave equations, path integral), there are three possible 
descriptions/definitions of the partition function: as a matrix element, as a solution to a system of 
linear differential equations (as an element of a -D-module), as a (functional) integral over trajectories in 
configuration and/or phase space (over field configurations). These very different definitions emphasize 
different properties of partition functions, and their equivalence implies deep and non-trivial relations 
and symmetries. Among these implications, there are integrability properties of partition functions, 
placing them into the class of generalized r -functions !2 a , which satisfy infinite sets of compatible non- 
linear difference-differential equations, (generalized Hirota equations). Moreover, partition functions 
of different models are often related by dualities and/or mirror-like symmetries. Despite their general 
definitions and symmetries, partition functions are rarely expressible through conventional special 
functions and possess complicated analytical properties, with all kinds of singularities and branchings. 

All this - universality, rich symmetry properties and impossibility to reduce the problem to known 
functions - means that the partition functions (r-functions) of string models are natural candidates 
for the role of the next generation of special functions: they should be studied without reference to any 
particular applications, their properties should be searched for and listed, particular interesting cases 
(at specific values of parameters, in asymptotics, on particular branches) should be enumerated, listed 
and tabulated - and finally collected in the special-functions reference-books. This task is even more 
natural, because, similarly to conventional special functions (of hyper geometric, elliptic and Riemann 
theta-functions families), all the r-functions are closely related to representation theory of Lie algebras 
and groups |2j. It is important to emphasize only that, while partition functions are usually some r- 
functions, inverse is not true: partition functions form a very special subclass in the set of r-functions - 
with additional structures and deeper properties (the situation is well modeled by the relation between 
Riemann and generic theta-functions: the former, representing truncated partition functions of free 
fields on Riemann surfaces, possess additional properties of crucial importance, e.g. the Fay identities, 
implied by the Wick theorem for the free fields). In other words, integrability theory is a natural part 
of the theory of partition functions, but only a small part, by no means exhaustive. 

In this paper we present a first iteration of the above-proposed analysis/listing for the simplest 
and very important string-theory special function: the partition function of the finite-size Hermitean 
1-matrix model (FSHMM) 0-0- The general scheme of reasoning is equally applicable to any other 
matrix model. The obvious next candidates to be analyzed from this perspective are quiver (conformal) 
matrix models [8]-|10| (see also the recent progress in the two-matrix (normal matrix) models 
and the Kontsevich r-function ^7j- After this is done, one would be capable to proceed to the 
geometric r-functions, associated with the topological sigma-models and models of quantum gravity 
in different backgrounds. However, of all partition functions, the matrix model ones are the natural 
candidates to begin with (see various aspects of matrix models in |18|-f20 j). Not only are they the 
simplest example, one can expect they are actually elementary building blocks for many other r- 
functions. An example of decomposition of a geometric partition function (the one for the CP n 
topological sigma-model) into a multilinear composition of the elementary building blocks (n + 1 
Kontsevich r-functions in this particular case) has been found by A.Givental |21|-|22| (see also [10 ). 
A simpler decomposition formula for the FSHMM (which, actually, appeared already in [22]) will be 
analyzed in the present paper. 
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• In the theory of conventional special functions one often distinguishes between the two levels of 
generality: one can look at generic solution to a differential equation and one can instead concentrate 
on particular solutions (branches), normally associated with specific integral representations (and/or 
particular representations of an underlying Lie algebra). To give a trivial example we can refer to 
the cylindric function (-D-module, associated with the operator z 2 d 2 + zd z + (z 2 — v 2 )) and its 
particular branches: Bessel and Neumann functions or Hankel functions. Any pair of these branches 
forms a basis in the space of solutions (in the -D-module) and can be fixed by choosing integration 
contours. 

The same is true for the string-theory special functions, particular (functional) integrals provide 
particular branches of generic partition function that is defined by a set of Ward identities (Shwinger- 
Dyson equations) - only the difference between branches can be somewhat more pronounced than in 
the elementary case. We shall elaborate on this in detail in the case of the FSHMM, beginning with 
a naive matrix integral in s|2j generalizing to arbitrary solution to the Virasoro constraints in ssOlBl 
(though today we do not have much to say when the genus expansion does not hold), and coming back 
to particular solutions: the Gaussian r- function (i) in s^and the CIV-DV r- function Znv[T\S](t) 
in sEJ With the help of the Givental-style decomposition formulas (s|2J), the CIV-DV r- function is 
constructed from a multilinear combination of the Gaussian ones, and, according to the arguments 
of s l4.Hl it can be considered as providing a basis in the linear space of all solutions to the Virasoro 
constraints (somewhat analogous to e ipx in the space of all the functions of x; in other words, just 
like e tpx ', Zdv[T\S] provides a kernel of an integral transform from functions of T into functions of 
S). Unfortunately, no feature is known yet which distinguishes this basis among other possibilities, 
although it has clear advantages of being provided by Givental-style decomposition formulas and being 
related to Whitham-Seiberg-Witten theory. 

• The goal of this paper is four-fold: 

- We formulate the task of tabulating the properties of matrix-model and geometric r-functions, 
defined as D-modules, i.e. as generic formal-series solutions to a system of Virasoro-like constraints, 
and make a preliminary list of those for the FSHMM. 

- We claim that description of the complicated phase (branch) structure of the r-function is a 
very important - and so far under-estimated and under-investigated - ingredient of this program, and 
make an attempt of such description for the FSHMM. 

- We introduce a hierarchical description of phases. First, with the help of the string coupling g, 
we introduce the gradation in powers of t- variables (coupling constants or times), and use it to select 
the phases possessing the genus expansion, where logarithm of the partition function (the prepotential) 
admits a semi-infinite expansion in powers of g. Second, a shift of finitely-many t- variables, t—*t — T, 
allows us to introduce different kinds of t-expansions (phases) with different quantities (actually, the 
differences between the roots of the polynomial of degree n, W'(z) = ^kTkZ k 1 = 
FliLiX- 2 — a i)) emerged in denominators. The polynomial W(<f>) of degree n + 1 plays role of the tree 
superpotential in the context of Dijkgraaf-Vafa theory |25j . Third, each of these phases is further 
split into infinite-parametric set with different positions of singularities (branching points) of the 
multidensities 1 

ml \ d \ 

pw(zi, ■.-,2 m )=n e fc+i n+ ^w(t) (i.i.i) 

where Zy/(t) = e^ 7 w is the FSHMM partition function. This family of phases is parameterized by 
an arbitrary function Fw(g) = F[g, T] = J-"w(t = 0) 9) of 9 an d n variables To, ... , T n _i. 

1 These quantities can be associated with the multiloop correlators (correlators of Wilson loops), /^FJ^Lj Tr ^'L^ y 

|26| (where the double brackets denote connected correlators) , when a matrix integral representation with some diagram 
perturbation expansion exists. The imaginary part of these correlators is naturally associated with multidensities of 
eigenvalue distributions in this case. This is why we usually call pw(zi, ■ ■ ■ , z m ) as multidensities. We also call them 
multiresolvent, since p(z) is a resolvent. 
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- We argue that an alternative approach to description of branches can be provided by the Givental- 
style decomposition formulas, representing a matrix- model/geometric T-function through an operator, 
acting on a product of "elementary r- functions" . In the case of the FSHMM, the "elementary r- 
functions" are the Gaussian ones, and the decomposition formulas give rise to the CIV-DV r- functions 
Zuv[T\S](t). Given W(z), the generic partition function, which is built from an arbitrary function 
Z[T] = e F ^ T ^ 9 with the help of the evolution operator Uw(t) (see s 14. 3(1 . 

Z w (t) = Uw(t)Z[T] (1.1.2) 
can be expressed through the CIV-DV r-functions, 

Z w (t) = J Z DV [T\S](t)C[S]dS (1.1.3) 

where the coefficients C[S] are defined by the expansion of "initial" Z[T] into a linear combination of 
Z DV [T\S] = e *W[T|S]/ fl * = z DV [T\S](t = 0), 

Z[T} = J Z DV [T\S]C[S]dS (1.1.4) 

2 Partition function as functional integral 

2.1 The notion of partition function 

• String theory studies partition functions 

S(<Po + <f>) 



ZM,S {t\4>o} = / # exp- 
JM 



M 9 

(1.2.1) 

S(0) = S o (0) + E%W) 
k 

(with complete set of vertex operators Vt{4>)) as functions of: 

(i) the space of fields Ai, parameterized by (f>, and equipped with some "natural" measure d<f>, 

(ii) the measure d(j,{(p) = d<j) e« S °^ on M. (composed from d<f>, the bare action Sq(4>) and the 
string coupling/Planck constant g), 

(iii) the infinitely many coupling constants and 

(iv) the background state (vacuum) fa. 

Usually the background field 4>q is a solution of the equations of motion 



6(f) 

but sometime their quantum counterpart, 

5dfi((j)) 



= (1.2.2) 

<t>0 



= (1.2.3) 

can be used, and/or a deformed action S{4>) substituted instead of Sq{4>) <|I-2.2|) and ()I.2.3|) (thus 
mixing the dependencies on </>o, g and tjA. 

The space A4, together with possible restrictions on the choice of measures dfi(4>), defines a par- 
ticular string model. 

The task of string theory is also to provide a unified construction of all interesting spaces Ai from 
some elementary building blocks (for example, from elementary simplexes). 

• Typically, 

M = Hom(W,T)/g (1.2.4) 
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is the space of mappings from the world sheet W into the target space T, factorized by some symmetry 
group Q. 



• For every A4 one can construct the new spaces A4g = M. ® Gr/U, where Gr is a group G in 
some representation R, factored out by action of its maximal subgroup U. The characteristic example 
is Gr = Mjv, the space of square N x N matrices (the iV-dimensional representation of GL{N) - the 
space of linear mappings of a vector space Vn into itself) , on which the unitary group U (N) acts by 
conjugations M — ► U~ l MU . This "branization" prescription (substituting numbers by matrices) is 
familiar from the studies of brane stacks [27] . 

For M of type (lL2~4l) . one can also switch to mappings between more general bundles, sheaves and 
complexes over W and T (substituting numbers by mappings of complexes, with complexes described 
by auxiliary variables like, for example, those in BW or BFW formalisms |28|-|33|). 

• For a given vacuum <po, it often makes sense to redefine the set of coupling constants in the action 
S((j)) in (|I.2.1|) : this is always possible if the vertex operators V^(ip) form a complete system (either 
strong or weak, see jSU for more details): 

5(0o + 0) = S (0o + </>) + E H V %^ + 0) = 

k (1.2.5) 
= 5o(0o + 0) + E4*° } ^ 

k 

The couplings tjf ^ are more adequate than for building up the perturbation theory around a 
particular background 4>q, but they themselves depend on (f>Q, and the Ward identities (like the Vira- 
soro constraints 
whenever tk = 0. 



soro constraints below) look more sophisticated if expressed through them. Of course, all 4 0o} = 



2.2 Matrix Model as the simplest String Model 

• The very starting example of string model is W = pt, T = pt (i.e. both the world sheet and the 
target space are just single points). The simplest associated brane stack provides M. = Mn/U(N), 
and the partition function is the one of Hermitean 1-matrix model of the eigenvalue type (Hj: 

z (™tr) it)= 1_ f eX p(l^t Tr^.....Tr^| (1.2.6) 

Vol 6r(iV) J Mn \ 9 ]? ) 

where the sum is over all non- negative integer- valued vectors k = (hi, . . . , k m ) of all possible lengths 
m. The measure 



N 



D<j> = J] #ij (1-2.7) 
is U{N) invariant, and, if cj) = U tpoU, where 4>d = diag(0j), 

N N 

Dcf> = DUl[d ( f> i l[( ( t> i - ( f> j ) 2 (1.2.8) 

i=l i<j 

DU is the Haar measure on Un, induced by the norm ||5C/|| 2 = Tr(U~ 1 5U) 2 , and the volume of the 
unitary group is 

. N-l 

Vol u{N) = DU ~ J] k\ (1.2.9) 

Ju(N) ^ 

The model is of eigenvalue type, because it is essentially reproduced by an iV-fold integral over the 
eigenvalues fa only. Integration contours for fa are not obligatory restricted to the real line (this can 
be considered as a sort of analytical continuation, see footnote |3J). 
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• Often instead of Zpf(t) the truncated partition function Z]y(t) is considered, with t? non-vanishing 
only for the vectors k of unit length, 



Z { ™ atr \t) 



Vol 



f Dc/> exp ( i V t k Tr fc I (1.2.10) 

u(N) Jm n yg fcgz+ J 

This is what we will use as a prototype of partition function of the finite-size (i.e. finite N) Hermitean 
matrix model (FSHMM). 

• A number of powerful methods exists for investigation of the truncated model ()I.2.10|) as an 
eigenvalue model [El El- The most important two are the method of orthogonal polynomials 3 a and 
the free-fermion formalism [51 E3- They can be used to derive the determinant representation [SUSIE] 
for Z/v, 

Z^ atr \t) = det< iJ<N n i+j (t) (1.2.11) 

where 

Hj = [ ^ exp ( - £ t k <\> k ) # (1.2.12) 
J \9 k=0 J 

In turn, it can be used to prove @] that Z^ is a Toda-chain r-function, i.e. satisfies an infinite system 
of bilinear differential equations like 



, f)2 7 {matr) ~ 7 {matr) r, 7 {matr) 
ryimatr) u ^ N u N N ryimatr-) ryVmatr) (r 9 1 q\ 

z - v —d% mi 8^ = Zn -i Zn ^ 1,1 



• In (jl.2.10j) the bare action is absent. As in every string model, it can be introduced by a shift of 
the coupling constants, t k — > t k — T k : 



^7 '(t) = TFT~ / d4>e~'" ww e- a ^ U9 (1.2.14) 

where VF(z) = J2k=i T k z k ■ Despite obtained by a shift of variables, this partition function essentially 
depends on W(z), especially on its power n + 1, since changing n completely changes structure of the 
perturbation series. In other words, the system undergoes a phase transition (the partition function 
is singular) at zero values of the coefficient of the leading (z n+1 ) term in W(z). 2 



3 Partition function as a D-module 

3.1 The notion of partition function 

• Partition function (r-function) is, in fact, a highly sophisticated quantity. Actually, it is not a 
function, but a formal D-module, i.e. the entire collection of power series (in t- variables), satisfying a 
system of consistent linear equations. Solution to the equations does not need to be unique, however, 
an appropriate analytical continuation in £-variables transforms one solution to another, and, on a 
large enough moduli space (of coupling constants t), the whole entity can be considered, at least, 
formally as a single object: this is what we call the partition function. Naively different solutions 
are interpreted as different branches of the partition function, associated with different phases of the 
theory. Further, solutions to the linear differential equations can be often represented as integrals (over 
spectral varieties), but integration "contours" remain unspecified: they can be generic chains with 
complex coefficients (in the case of integer coefficients this is often described in terms of monodromies, 
but in the case of partition functions the coefficients are not restricted to be integer). A model of 
partition function is an integral formula which has enough many free parameters to represent the 
generic solution of the differential equations in question. 3 

2 This situation is much similar to matrix models of the Kontsevich type, which also drastically depends on the matrix 
model potential, especially on its power |16|. 

3 Several familiar examples can be helpful to clarify these notions. 
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3.2 Virasoro constraints for the FSHMM 



• Coming back to the FSHMM, its original definition has been done in terms of a matrix integral 
()I.2.10|) . This integral depends on (infinitely many) parameters tk and, according to above arguments, 
implies another definition of the partition function: as solution to the (infinite) system of consistent 
linear (differential) equations (Ward identities or Picard- Fucks equations), i.e. as a D-module. In the 
particular case of FSHMM, these equations, 

L m Z(t) = 0, m > -1 (1.3.10) 



d , ^ d 2 



Lm ~^ ktk 3^Z + ? ? dtM (1.3.11) 



k>0 UZ k+rn a+b=m 
a,b>0 



First, a flat SU(2) Seiberg-Witten modulus a(A) |35| satisfies 

(4A(A - l)dl + l) a(X) = (1.3.1) 

possesses integral representation, 

a(X) = f J <fa (1-3.2) 



j c V x(x-l) 
and the model is provided by the generic contour 

G = a[0, 1] + /3[1, A] (1.3.3) 

which is the generic linear combination (with complex(!) coefficients) of A and B cycles on the toric spectral curve. 
Alternative description of the same object (the toy partition function) is in terms of the Seiberg-Witten 1-form 



dS(x,\) = J j \ dx (1.3.4) 

y x(x - 1) 

related to the original equation 1)1.3.1^ through 

(4A(A - 1)^ + 1) dS(x, A) = d\J x{ *~^ S (1-3.5) 

(where = means cohomologically equivalent). 

Another example already discussed in the introduction is the cylindric functions. Their defining equation is 

[\ 2 8l + \d x + (A 2 - v 2 )] Z V (X) = (1.3.6) 

and an integral representation is 

Z J X ) = 1- [ e -i-Bine+iu0 de (I 3 7) 

27V Jc 

The model is given by the generic linear combination of two contours, say, chosen as in 8.423 of |24] (this choice fixes as 
the basis the Hankel functions). 

One more example is the Airy function (a simple example of Kontsevich model |12 | -|17 | ~). with the defining equation 

(d\ + A) A(X) = (1.3.8) 

and an integral representation 

A(X) = J exp + A^ dx (1.3.9) 

This D-module is a particular case of the cylindric functions with v = i. 

4 The origin of these equations is invariance of integral 11.2.101 under the change of integration variable <j> — » (f>+e m (f> m+1 
|36j . Similarly, integral ljl,2.6fl is invariant under the changes (f> — » <fr + e m g<^) m+1 Tr cj> kl • . . . • Tr 4> k ' . The resulting set of 
constraints is a corollary of 11.3.1011 and the obvious relations 

dZ _ d m Z 
dt% dt kl . . . dt km 

The inverse way, from the equations 11.3.1011 to the integral representation 11.2. 10L is discussed in [Sj. 
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are called |37j Virasoro constraints, since L m ,L n = {m — n)L m+n (in fact, this is only the Borel part 
of the Virasoro algebra). 

From now on we define the FSHMM partition function as the generic solution to the system 
()I.3.10|) . so that the matrix integral ()I.2.10|) is its particular integral representation. One of our tasks 
below is to analyze how general this integral representation is and whether it can be treated as a model 
of the FSHMM partition function. 



Since 



9 L 



r (1.3.12) 

dt 

one can diagonalize the operator d/dto and impose an additional constraint on the partition function 

= SZ S (1.3.13) 

ot 

with constant S. The integral representations ()I.2.10|) and (jl.2. 14j) obviously satisfy this relation with 
the integer S = N. The generic solution to the Virasoro constraints can be decomposed into a linear 
combination of Zg's with different constant S's, but, in practice, this decomposition does not look 
very natural. 5 

Note that, for particular values of N, the partition function Z^p atT ^ satisfies additional differential 
equations, impled by the fact that Tr</> fc with k > N can be expressed through an algebraic combination 
of Ti(j) k with k < N. For example, for N = 1, 



d_ 



Z { ™{ r) = (1.3.14) 



We do not impose these constraints in our definition of the generic FSHMM partition function: the 
only requirements are (|I.3.10|) (and (jl.3. 16|) . to be defined below). 

• Solutions to the system ()I.3.10|) can be looked for among the formal series in powers of t^, and 
these series are different, depending on what is allowed to appear in denominators. Different such 
possibilities are interpreted as different phases of the matrix model or as different branches of its 
partition r-function Z(t). Formally, they can be regarded as different approximations and analytical 
continuations of the integral (|I.2.10|) . These different phases are labeled, in particular, by the bare 
actions, i.e. by consideration of (|I.2.14|) instead of (|I.2.10|) . 

• Introduction of the bare action W(<fi) = J2k=i Tk4> k modifies the Virasoro constraints (|I.3.10|) in 
a trivial way, by shifting 

d 

L m ^ LZ = L m -J2 kT kJf (1.3.15) 

fc>0 ° l k+m 

but the solutions Z\y(t) = e^^^ 2 can be now looked for among the formal series in positive powers 
of the t-variables with T's emerging in denominators. Also, now the set of the Virasoro equations 
should be supplemented by additional constraints 

^ + f^ = 0, Vfc = 0,...,n + 1 (1.3.16) 
oT k dt k 



5 What is worse, while the constraints Ni = const in the Givental-style decomposition formulas sect. II. 2 are clearly 
generalizations of 11.3.131 (and, moreover, J^iLi — N = S), no operators are known that have Ni as their eigenvalues. 
This causes problems with justification of linear decomposition of every branch of partition function into the Givental 
functions 111.2.511 with constant iVVs. 
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3.3 Virasoro constraints and loop operators 

• A powerful technique for dealing with the Virasoro equations (jl.3. lOf) is provided by the loop operator 



V W = E 7FTT 

fe>0 



z k +! dt k 



(1.3.17) 



and the two complementary projectors PJ onto negative and non-negative powers of z, 





+oo 






E 


= E t k z k 




k = — OD 


k<0 


Pi 


+oo 

E ^ k 


= E ^ fc 




k=—oo 


fc>0 


With the help of these operators and the 


generating 


function 




oc 





(1.3.18) 



v{z) = J2 M* 

k>0 



(1.3.19) 



one can write down the entire set of the Virasoro constraints 1)1.3.10(1 in the form of a single loop 
equation 

P- W(z)p(z)] + p\z) + g 2 V(z)p(z) = (1.3.20) 

for the resolvent (density) 6 



1 BT 



k>0 



z k + l dt k 



Ea. ()I.3.20|) should be complemented with the zero-curvature condition 

V(x)p(z) = V(z)p(x) 



(1.3.21) 



(1.3.22) 



After the shift ()1. 3.15(1 . the modified Virasoro relations for 

pw{z) = V(z)J : w 



arc 



(1.3.23) 
(1.3.24) 



u(z), 



P- [W\z)p w {z)\ = p 2 v (z)+g 2 V(z)p w (z)+p- [v'{z)p w {z)] 

In considerations below we make use of one more loop operator, defined for arbitrary function 

u'{z) 



R u (z)=J2^ 

l>0 



z 



l+l 



d_ 

dti 



(1.3.25) 



6 The Virasoro constraints 11.3.101 can be also realized in free field terms. Indeed, consider the {/(l)-current d z 4>{z) in 
the theory of one free scalar field realized in terms of Heisenberg algebra acting on the space of functions of time- variables 

ifc> 

d^{z) = ^v'(z)+gV(z) 

Then, the Virasoro generators are given by the standard expansion of the energy-momentum tensor 

Now the loop equation is just the requirement that the partition function as a function of time variables realizes the 
vacuum w.r.t. the Virasoro algebra: 

P~T{z)Z = 

This language proved to be very effective in dealing with different quiver models and with continuum limits, see |38l 15). 
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Characteristic property of this operator is that, for arbitrary function g(t) of time- variables, we have 



u'{z)V{z)g{t) = R u {z)g{t) 



(1.3.26) 



In particular, if u(z) is the potentials v(z) or W{z), and g(t) is F(t) = g 2 log Z or F w (t) = F(T, t) 
g 2 log Z\y , one obtains for the corresponding pieces in the loop equations 1)1.3. 20JI and ()I.3.24|> , 



P- [v'{z)p{z)\ = v'{z)p{z) - P+ [v'(z)p{z)\ = v\z)p{z) - R v (z)F 



(1.3.27) 



and 



P- [W'(z)p w (z)} = W'{z) Pw {z) - P+ [W'{z) Pw {z)] = W'(z) Pw {z) - R w {z)T w (1.3.28) 



For polynomial W(z), the last term at the r.h.s. is just a polynomial of degree n — 1 in z. Moreover, 
due to (11.3.16)1 . one can also write 



R w {z)Fw = ~ ^2{k + l + 2)T k+i+2 z 

k,l>0 



dT, 



(1.3.29) 



Hereafter, we use ~ to denote operators containing T-derivatives, while those containing i-derivatives 
are labeled by * . 

With the help of R\v(z), we can rewrite p. 3. 24)1 in the most convenient form 



where 

Note also that 



W'(z)p w (z) = p 2 w (z) + f w (z) + g 2 V(z) Pw (z) + P- [v'(z) Pw (z)] 
f w (z) = P+ [W'{z)pw{z)] = R w (z) log Z w 



V{x)v{z) = P- 



z — X 



V{x)v'(z) = -P- 



1 



and for any function h(z) such that P^~[h(z)] = 0, 



V(x)P" [v'(z)h(z)] = -P- 
where d x = d/dx. The commutator 



h{z) 



V(x),R v (z) =d x P, 



(z — x) 2 
V(z) 







(z — x) 2 . 
h(z) — h(x) 



z — x 



z — X 



(1.3.30) 
(1.3.31) 
(1.3.32) 

(1.3.33) 
(1.3.34) 



3.4 Solving Virasoro constraints recurrently 

• In this section we assume that Zw(t) possesses genus expansion, i.e. that 

log Z W (t ) = g- 2 T w (t) = £ g 2p - 2 f$ (t) 

p>0 

(see s !4.1l for discussion of this assumption). 

Actually, the partition functions can be obtained by one of the two procedures: 



(1.3.35) 



- Virasoro- constraints method, i.e. by solving the loop equations 1)1.3.30)1 . and 

- Givental method, i.e. by making use of decomposition formulas ()II.2.5)1 . to be introduced below 
in s.110 

In the framework of the Givental method, an arbitrary phase of the FSHMM with given W(cp) 
of degree n + 1 is rather described in terms of n Gaussian models, according to the decomposition 
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formula (111.2.51) . This means that the averages depend, in fact, on the additional set of variables iVj, 
characterizing a particular saddle-point (vacuum) of the FSHMM. We go into more details along this 
line in s.II|2] below, while here we concentrate on analysis of the Virasoro constraints ()I.3.30|) . which 
provides an a priori richer variety of possible solutions/phases (though they would still have to be 
linear combinations of the Givental-style ones). 

• In order to convert (|I,3.30|) into a solvable set of recurrent relations, we expand 

p w (z\t) = g 2 V(z) log Z w (t) = V{z)F w {t) (1.3.36) 

in powers of g 2 and t's 

pw(z\t) = ^2 ^ j---jv{zi)...v(z m )p i $ m+l \z,zi,...,z m ), 

p,m>0 

f W (z\t) = ~ f ■ ■ ■ f v ( z l) ■ ■ ■ v(z m )fw m+1 \z\z 1 , . . . ,Zm) 



(1.3.37) 



p 2p 

p,m>0 

In this way, we introduce the full set of multidensities 

p$ m \z 1 ,...,z m ) = V(z 1 )...V(z m )F$\ t=o (1.3.38) 

and auxiliary polynomials (which distinguishes between different phases for a given W{z)) 

f$ m+1 \z\ Zl , ...,z m ) = R w (z)p$ m \z u . ..,z m ) (1.3.39) 

In fact, ea. (jl.3.39|) guarantees that the constraints (|I.3,16|) are fulfilled. 

Acting on ea, H1.3.30|) with the operator V(-?i) . . . V(z m ), making use of (|1.3.33|) and putting all 

tk = afterwards, we obtain a double-recurrent (in p and m) relation for the multidensities p^ 

W'(z)p$ m+1 \z, z h ...,z m )- f$ m+1 \z\ Zl , ...,z m ) = 

= ^ ^ Pw \ z : ' Zi i ' • • • ' ^imj 'Pw \ z i z jn ■ ■ ■ ■> z jm 2 ) + 

g m\+rri2=m 



m a (p\m) / _ \ (p\m) / \ 

_l_ ° Pw \ z i z li ■ ■ ■ i z ii ■ ■ ■ i z m) — P\\r \Z\,...,Z m ) 

■ ■x dzi z — Zj 



(1.3.40) 



+V(z)p w -^ m+1 \z,z 1 ,...,z m ) 

Together with (|I.3.39|) this relation is enough to find explicit expressions for all the multidensities 
through W(z) and f w ^\z). In fact, the latter polynomials (all of degree n — 1) are not independent, 
since (|L3~39l) for m = 0, 

f^ 1 \z) = R w (z)F^[T] (1.3.41) 

holds only because of the zero-curvature constraint (|I.3.22|) and expresses all the /-polynomials through 
a single function of T (i.e. of W{z)) and g, the prepotential at t = 0, 

oo 

F w (t = 0,g)= F[g,T] =Y / g 2p F^[T] (1.3.42) 

p=0 

• The T n and T n+ \ dependencies of F[g, T] are not arbitrary, since the and L w Virasoro 
constraints can be consistently truncated to t = and then allow one to express dF / dT n+ i and 
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dF/dT n through dF/dTi with I = 0, . . . , n - 1, 



dF (0) 1 | ^ dF (0) 

dT n ~ (n+l)T n+1 ^ [L+ 1)11+1 dT, 



2 



(1.3.43) 



dT n+1 ~ (n + l)T n+1 1, dTo J 
~(n + l)T„ + i g /Ti ^T " (n + 1) 2 T 2 +1 + 

• As an immediate corollary of (|I.3,40|) . we obtain for p = and m = 

p (Ji) (z) = W^LM!) (1.3.44) 

with 

y 2 / (z) = (W'(z)) 2 -4/g |1) (z) (1.3.45) 



Sometime it is useful to represent the polynomial fyy^iz) through peculiar n parameters Si |39j 



4/fW = ^)E^T (I-3-46) 

1=1 1 



where cti are the roots of W'(z) = JJ^ =1 (z — cti) (so that there are no poles at the r.h.s. of (|I.3.46|) ). 
If S's do not depend on T's, the corresponding truncated prepotential is especially simple, 

n 

F^[T\S] = -Y / S i W(a i ) (1.3.47) 
Indeed, since W'ipn) = 0, the second term at the r.h.s. of 

1=1 8=1 



M WiE^, A/; = W"(ai) (1.3.48) 



vanishes, and the action of Rw(z) straightforwardly converts the first term into the r.h.s. of (|1.3.46|) . 
See Tables below for more explicit expressions for the multidensities. 



4 Comments on the number of solutions 



In this section we discuss an a priori way to introduce phases, and naturality of their labeling intro- 
duced in the previous section. We also emphasize that an interesting task is to construct a matrix- 
integral model of the partition function defined by the Virasoro constraints (|I.3.10|) or (jl,3.20|) . and 
briefly discuss in what sense the Givental-style decomposition formulas (which are related to matrix 
integrals and give rise to the CIV-DV partition functions) do provide a solution to this problem. 

As we already saw above, an important ingredient 7 in the theory of matrix model integrals, is 
relevance of the shifts of t-variables: Z(t — T) and Z(T' + 1) with different T and T' can be different 
branches of the same partition function, moreover, they are further ramified into (generically infinitely 
many) sub-branches. Below we explain the meaning of this claim and describe the pattern of branches 

7 This is new as compared to the simple example of 11.3.11 . but can be already observed in the integral representation 
for the cylindric functions 11.3.71 . Indeed, the integration contour, say, for the Hankel functions is quite fancy so that 
the shift of integration variable generates from any one Hankel function their linear combination. 
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and their properties. 8 The first question to address is to enumerate/label the branches, i.e. different 
formal series solutions to the system 3.10|l . As already mentioned, an origin of complications is that 
the series depends on ratios of t-variables and different branches occur when different combinations 
of t's are allowed to appear in the denominators. One can distinguish three levels in classification of 
branches, characterized by the parameter g, W(z) = J^TkZ k and fyy (zi, . . . ,z m \g,T). 

4.1 Types of g dependencies 

• The ^-dependence characterizes the degree of singularity when all t/. — > or all — > oo. From the 
very beginning one can distinguish between three possibilities: 

- Expansion of Z in non-negative powers of g. The leading term of expansion in this phase describes 
the "naive g = limit" and satisfies the equations 



1^ kt kTr, = 0, m > -1 (1.4.1) 

fe>o dtk + m 



- Expansion of Z in non-positive powers of g. In this phase the leading term of expansion corre- 
sponds to the "naive g = oo limit" and satisfies the equations 

a,b>0 " " 

a+b=m 

— Alternative to these two cases is a Laurent expansion in powers of g, infinite in both positive and 
negative powers of g (if the set of powers is semi-infinite, multiplication of Z by appropriate power of 
g puts it into one of the two above classes). Infinite Laurent series can be sometimes further separated 
into classes, if, instead of Z, one considers the prepotential J- = g 2 log Z . Then 

= £ 9 2P F$ (1.4.3) 

q<p<s 

gives rise to an infinite Laurent series for Z whenever q < and s > 2. Different pairs (q, s) give 
rise to different branches of Z. In fact, solutions to (|I.3.10|) exist only if q = — oo or s = oo. If both 
q = — oo and s = oo, one can repeat the trick and take, say, log T instead of J- and split the "point" 
(—00,00) further and so on. Also p in (|I.4,3|) can be, at least, half-integer (more general expansions 
in non-integer powers of tk, though quite interesting, are even more exotic). Of all allowed T q , s -, the 
most attention in the literature is paid to T = J-0,001 usually with integer (rather than half-integer) 
p. In this case (|I.4.3|) is called genus expansion. J-^ is the spherical contribution to the prepotential 
T = ^0,00) P labels the genus. The half-integer p's, if included, are associated with Riemann surfaces 
with boundaries or with non-orientable surfaces (both have doubles, which are ordinary Riemann 
surfaces with integer genus). 

• In the most of phases J- qs , for given q and s (in particular, ^~o,oo 

) possesses additional free pa- 
rameters, which are not fixed by the Virasoro constraints (|I.3.10|) . In s.3.4 we labeled these parameters 
by / (actually by f(z) made from an arbitrary function of T- variables F[g,T], (|I.3.41|) ). Every linear 
combination 

Z {°} = E c q,s;f ex PFq,s;f(g\t) (1.4.4) 

with ^independent (but, probably, ^-dependent) coefficients c qtS -j is still a solution to the linear system 
(|I.3.10|) . and c q:S -j can not be absorbed into redefinition of /. These coefficients c q:S -j are analogues 

8 After this is done, one can switch to studies of the continuum limits of these branches, which arise when N — > 00 and 
other parameters behave somehow, thus providing (depending on behaviour of these infinitely many parameters) 00 °° 
many types of continuum limits (there can be rather few on the lowest branches - and these are all so far described in the 
literature). However, below we stop at the level of finite iV: classification of branches is a precondition for classification 
of continuum limits, and it is already enough sophisticated. 
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of a and (3 in (|I.3,3|) and can be interpreted as substitution of the integration contours by cycles with 
complex coefficients. In what follows we pay no more attention to the freedom associated with c q ^ s -j 
and concentrate on the D-module generators Z q)S -j = exp Tq iS -j{g\t)- Our next task is to explain the 
origin of the free parameters f. 9 



4.2 The shifts of t and the origin of /-parameters: an oversimplified example 

• Once the ^-dependence is specified, one can start enumerating solutions to the system (|1.3.10j) . To 
describe different phases, we shift the variables t k — > —T k + t k and assume that only the first n + 1 of 
Tfc's can be non- vanishing. 

In order to understand how the new free parameters appear for non-trivial W(z), it is instructive 
to start from the "naive g = limit" , i.e. from the partition function satisfying eq. (|I.4.1|) . Note that 
it has nothing to do with the partition function in "spherical" (genus p = 0) approximation, exp 
Instead of (ll.3.2fi|) . we have for 

p (z)=V(z)logZ (1.4.5) 
(note the absence of the factor g 2 in this definition) evaluated at t k = (i.e. when v(z) = 0) 

P- [W'(z)po(z)] = 0, i.e. p (z) = J£jL (1.4.6) 

where deg{f{z)) < degW'(z), since the operator V(z) and po(z) in (|1.4,5|) are expanded in negative 
powers of z only. This means that f(z) is an arbitrary polynomial of degree n — 1. In other words, 
we see that, for given VK(z) of degree n + 1, we get an n-parametric set of solutions, even if all tk = 0. 
Switching on t^s is equivalent to raising the degree of W(z), and thus could lead to a new arbitrariness, 
to new arbitrary polynomials like f(z). However, the coefficients of these polynomials are actually 
functions of T- variables restricted by additional constraints (|I.3.22|) . The formalism, developed in s !3.4l 
below, allows one to demonstrate that all the /-functions are made from a single arbitrary function 
F(T) of n variables To, . . . , T n _i . 

In our oversimplified example this can, of course, be also seen without any special formalism, just 
by looking at the truncated Virasoro constraints 

g Mi ^W = f;« t ^W, ra >_i (i.4.7) 

fc=l UL k+m fe=1 UL k+m 

Putting all tk = 0, we express all the derivatives n with k > n through those with k < n. Now 

take the i^-derivative of (|I.4.7|) . 



*=o 



^ d^ w (t) MV(g) " dT w {t) 

l^ kT k7T + Tu =l ~tu V 2^, kt k^l 5T> m>-l (1.4.8) 

k=1 ot k+rn oti oti +m k=1 ot k+m oti 

Putting all t k = 0, we express all the second derivatives with k > n through those with 

k < n. Using the symmetry between k and I we can conclude that, of all the second-derivatives at 
t = 0, only »,'^f with k,l < n are left unrestricted by the Virasoro constraints. Continuing 
in the same fashion, one proves that the truncated Virasoro constraints fix everything except for 



dt kl ...dt kp 

truncatec. 



with k\,...,kp<n, with any p. This means that, indeed, the only freedom left after the 

t=q 

Virasoro constraints are imposed is in the choice of arbitrary function F[T] = J~w(t = 0) 
of n variables Ti, . . . , T n -\ (since the dependence on T n and T n+ \ is fixed by the constraints L^it = 
0)Z w (t = 0) = L^(t = 0)Z w (t = 0) = 0, see ICT ). 

9 The separation of "degrees of freedom" into / and c/ is not absolutely obvious. We described a possible approach to 
the separation with the help of basis in the Hilbert space of solutions in the very end of introduction, and more details 
about the corresponding evolution operator are given in s 14. 31 below. 
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• Just the same is going to happen when g ^ with more interesting prepotentials Tyy . They 
depend on the choice of W(z) and on additional functions like f(z). All these functions, which 
are getting more and more complicated with increasing genus p and with increasing complexity of 
correlators (what is equivalent to increasing the number of non- vanishing, but infinitesimal variables 
t k ), are, in fact, made from the T-derivatives of just a single arbitrary function of To, . . . , T n _i, that 
is, the truncated (evaluated at t k = 0) prepotential F(T\g). 



4.3 Evolution operator and basis in the Hilbert space of solutions 

In this section we claim that, for any W(z), i.e. for any choice of T- variables, there is an operator 
Uw(t), which converts any function Z[T] of To, ...,T n _i (with prescribed dependence on T n and 
T n+ \) into Zw(t) = Uw(t)Z[T], which satisfies the Virasoro constraints, L^(t)Zw(t) = 0, m > 
— 1. Moreover, if Z[T] = J2a°aZ^[T], then Zy/(t) = J2a°aZ^ (t). The construction is, in fact, a 
straightforward generalization of that in s 14.21 to the case of g ^ 0. 

Namely, for given T's , we make use of the Virasoro constraints L^Z(t) = 0, i.e. 



n+l o y 

E kT k - ° Z 



, , dz \- 

E^^t— + E 



d 2 z 



dtk+m k=1 dt 
and their multiple i-derivatives to recurrently express 



k=1 —k+m k=1 vt k+m a+b=m dt a dt b 

dPZ 



m > — 1 



(1.4.9) 



dt kl ..M k% 



with all < ki < oo through 



J2s d p+s z 



p+s 



t=0 



t=o 



with all < lj < n 



d p Z 



dt kl . . . dt kp 



EJ2sT>h—h+' 
9 u k 1 ...k p 



t=o 



Q P +S Z 



dTu ...&Ti 



p+s 



t=o 



(1.4.10) 



E 



Js-n'i-'p+s 
9 u k 



t=o 

contains a , a , a . 

f=0 dt a dt b dt c 



t=Q 



4=0 ' 

this time with 



This is a little more tedious procedure than in the case of g = 0, but still it is straightforward, and the 
sum over s is finite, from to the integer part of the expression for contains QiQ th 

but with a, b < k — n — 1, further, the expression for ^ g t 
a, b, c < k — 2n — 2 and so on. 

Now we can define the operators 

k\ . . . 



t=0 



(1.4.11) 



kl - k ? dTu...dTi. 

s b l i p+s 

0<h,. ..,l p + s <n—l 

Because of their origin in the Virasoro constraints, these operators satisfy some obvious relations. 
Actually these are the eas, (|I.4.14j) below. Now we construct the evolution operator Uw(t) as a series 
in these T>-operators 

1 - (2) 1 
fktiD™ + —i 

The fact that for any Z[T] 



u w (t) = i + t k b { k 1] + U k tti)$ + ^WmD^ + ... 



L%(t)Z w (t) = L%(t)U w (t)Z[T] = 



(1.4.12) 
(1.4.13) 



or, simply, that 



i%{t)u w {t) = ( e kT k b^ m + £ b% I + E *« I <l + E ^gUi + E ] + 



n( 2 ) 



r)(l) 



n( 2 ) 



,(3) 



a+b=m 



a+b=m 



n( 2 ) 



6(3) 



a+b=m 



(1.4.14) 
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is equivalent to vanishing of all the linear combinations of operators in brackets, and these are the 
characteristic equations for the .D-operators. Note that we proved that the evolution operator is the 
same for any values of the arbitrary parameters / (or for any function Z[T]) once W(.z) is fixed. 
This proves that if Z[T] = T,a c aZ {a) [T], then Z w (t) = T, a c a Z$(t), and means that "orbits" of the 
evolution operators are completely parameterized by W(z). 
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Part II 

Particular Special Functions 

1 Gaussian partition function 

• The Gaussian partition function (phase) is distinguished by existence of diagrammatic representation 
for correlators and thus by the possibility of interpreting all the formulas in terms of combinatorics of 
fat graphs. This makes the Gaussian partition function an important special function by itself, even 
without a reference to the generic FSHMM partition function, of which it is a particular branch. More- 
over, the Givental-style decomposition formulas (see s.IlEJ) represent various branches as polilinear 
combinations of the Gaussian partition functions. For these reasons, the Gaussian partition function 
Zq (t\N) deserves a special attention, and we consider it separately, though a significant part of the 
formalism below is a particular case of the generic formalism for non-Gaussian models. 

• The Gaussian matrix model (phase) is the one with W(x) = yj; 2 (i.e. T k = ^-5 kt 2, where M is 
a number), and the saddle point (see s.I.2.1) 0o = 0, 

z£f(t\N) = 77^— / ^exp f-^TrcA 2 + - £ t k TrA = 
Vol u(N) J \ 2g g^ Q J 

2vrs\ y t kl -..t km / kl k v 

M J 4^ g m m\ \ I G 

k 

Here the brackets denote averaging over the random matrices with the Gaussian measure, exp ^— ||Tr^ 2 
Each correlator at the r.h.s. is a sum over the Feynman-t'Hooft diagrams {fat graphs) with m 
vertices of valencies k±, . . . , k m . Each vertex contributes the factor of g , each propagator, denoted 
by a double-line (strip of the fat graph), contributes the factor of g/M, and each closed line (boundary 
component of the fat graph) contributes the factor of N = S/g, where S = gN is the t'Hooft's coupling 
constant 10 . The prepotential Tq = 9 2 log Zq is a similar sum over the connected diagrams only. If 
we switch from the N, g variables to S, g, then particular terms of the expansion 

oo 

r G {t\S) = Y J 9 2p ^\t\S) (II.1.2) 
are sums over connected genus-p fat graphs 

(p) 

(i! -^ir r ....ir r; | = 

(II.1.3) 



J%\t\S) = F<?\S) + ^ ' ' ] tkm ((^ • • • • Tr ^ 

k 

= Fg\s) + E ^ / ■ ■ ■/«(*!) • • • v(z m )p^ m \ Zl , ...,z m ) 



m>l 



10 We remind that it can make sense to assume that S can still depend on g 

oo 

S = S k g k 

k = 

However, the expansion in (111.1.21 is defined for S = gN, otherwise one could write 

oo 

^a(t\S) = J2 92P:P o(t) 



p=0 

with 

,: >"-./,; 



(t) = (t\S ) + Si^^ 0) (t\S ) 



etc. 
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where the i-independent quasiclassical prepotential 

= - y log ^6 P , - (log VoLy^) 00 (1 - 6 Pi0 ) (H.1.4) 

and the Gaussian multidensities are 

p^ |m) (z!, ...,z m ) = {{T{cj>\z x ) . . . T{cj>\z m )))f (IL1.5) 

with T(<f>\z) = Tr-j^r, and ((. . ■))^ p ' denoting contribution of the connected Feynman-t'Hooft diagrams 
of genus p to the Gaussian average with the g and N -independent weight exp (— 4rtr0 2 ^. Formally, 
the size of the matrices is now S, not N: tr I = S (accordingly, we change the notation Tr for tr). 
This calculus is justified by the Euler's theorem stating that 

— # of vertices # of edges —# of loops 2xgcnus— 2 Qj i g\ 



• The Gaussian multidensities pQ^ m ' can be obtained recurrently in m and then in p by solving 
the Virasoro constraints (|I.3.10|) rewritten in terms of the prepotentials T 



M w 3 ~= E E E ^r^^m 

a + b=m<j+r=p " a + 6=m " fc>0 K-fm 

a,fe>0 <?,r>0 a,6>0 

m > —1 

Summing over m with the weight z~ m ~ 2 , one obtains for the generating function 



_1 ftF G (tH 
fc>0 r 



E^ E ~\ f ■■■ jv(z 1 )...v(z m )p { ^ m+1 \z,zi,...,z m ) 



(II.1.7) 



(ILL 



p=0 m>0 

the equation (c/. (|I.3.3Q|) ^ 

M(zp G (*|t) - 5) = p G (z|i) + g 2 V(z) PG (z\t) + P- [t/(z)p G (s|t)] (ILL9) 

where 

P- [v'(z)p G (z\t)] = v'{z)p G {z\t) - {{^ V '^\~_ V ' z {Z) )) g (n.1.10) 

while 

P- [zp G {z\t)} = zp G (z\t) - (tr /) = zp G (z\t) - S (ILl.ll) 

Putting v(z) = in equation (|II,1.9|) . one obtains a recurrent relation in p for the densities. Similar 
recurrent relations for multidensities with m > 1 are obtained by acting with V(zi) . . . V(z m _i) on 
flll.l.9JI . making use of (jl.3.33|) and putting v{z) = afterwards. This double recurrent procedure 
(with recursions in p and m) allows one to obtain every particular p( p \ m \ See Tables below for some 
explicit expressions. 

• The first (genus-zero, single-point) density, satisfying 

M(zp^ 1) (z)-S) = (p ( S ll) (z)) 2 (IL1.12) 
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is the celebrated semicircular distribution 18 , 



P ^ (z) - — I z - \ z 2 -—) - MZ - VG{Z) - T ^ - - + — + (II 1 13) 

It is a generating function for the Catalan numbers 

4*T(fc + l/2) 

Cfc = r(i/2)(fc + i)! (ILL14) 

counting the numbers of planar leaf diagrams with k leaves. 
The two-point density satisfies 

(Mz - 2 P ^\z)) P ^\z,x) = d /o^-p^Hx) (H.1.15) 

so that 

(0|2), 1 ( - 4M5 



2(z - x) 2 V yG{z)yG{x) 

J_ (-1 + J££_^ + 1 Mz 2 - 2MS - zy G {z) + Q / _L_ 
2x 2 V Vg{z)) x 3 y G (z) \x 2 



(III. 16) 



The first two terms at the r.h.s. are the generating functions of the correlators 

^ «tr* tr/))™ and «tr^ t^))™ (II.1.17) 

The latter generating function can be alternatively obtained by taking the derivative of p G ^ over 
T2 = M/2 provided that 

dS 

—— = for the Gaussian model (II. 1.18) 

0T2 

(note that this is not true for the generic FSHMM of s. 1 13. 41 even for the quadratic potential W(z) 
2 

• Along with the Gaussian multidensities 



^i—: whenever S 7^ const the diagram technique breaks down even for the quadratic potentials). 



*»> = ((^-^)f <IIX19> 



G 



one can consider their inverse Laplace transforms 



r$ m \x u r,„ ) f... f e x ^+- + ^ P ^ m \z 1 ,. . . , z m )d Zl ...dz m = ((Tt . . . Tr e*"*))^ 

(II.1.20) 

and their generating functions like 



-,2p f f poo poo 

dx\ . . . dx m x 



e xz p(z\t)dz = } — - <$>... <$> dzdz\ . . . dz m / . . . / 

Pi tlo m!/ 1 Jo Jo (ILL21) 

xe - XlZl -...- XmZm - {zi) jj^Kl)^ XI,..., x m ) 
The direct Laplace transform converts rj G (x) back into p G (z) 

POO 

PG {z\t) = / e- xz VG (x\t)dx (II.1.22) 
Jo 
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and 



pg m \z 1 ,...,z m ) = / .../ e-^-- x ^ V ^ m \x 1 ,...,x m )dx 1 ...dx m (II.1.23) 







• The (one-point) Gaussian densities satisfy the sum rule, calculating the total number of "flower" 
diagrams contributing to (Tr(/) 2fc ) 

pj z 2k p (p\l) {z) = {2k _ 1)H ^ (ILL24) 

The sum rule 1)11.1 .24(1 is a particular case (at N = 1) of a more general identity |4U1 141 j . providing 
a generating function for all (any genus) Gaussian one-point functions (densities) 

,r , T 2 2 (2N 3 + N)X 4 
AT _|_ J\T X -I- h 

^ (2k - IV ""^ /u 2r 2 I \ 1 - r 2 / * / 3 
(5iV 4 + WN 2 )x 6 + (14iV 5 + 70iV 3 + 21N)x 8 + (427V 6 + 4207V 4 + 483/V 2 )x 10 ^ 




15 105 945 

(fc+l/2 

r(i/2) 



Here (2fc - 1)!! = 2 k r §tj£ ) so that (-1)!! = 1. 



• The Gaussian partition function Zq(J,\N) can be also represented as a member of the family of 
Kontsevich models |12j-jl7j. [H]: it coincides with the Gaussian Kontsevich model with logarithmic 
potential 

J dHNxN6 = (27F) fdX nxn e-K^ (ILL26) 

where the time variables t^s are expressed through the external matrix in the Kontsevich integral, 
ifc = ^TtH~ k + ^5fc 5 2) the size of matrix in the r.h.s. is n x n and in the l.h.s. is N x N, and 
these parameters are absolutely independent. The multidensities p( p l m ) can be also reproduced from 
quadratic Kontsevich model (we postpone their discussion for a separate publication devoted to the 
Kontsevich type matrix models). 

• The genus-zero three-point functions, in particular p^ 3 \z\, 2 2 , 23), are expected to satisfy the 
WDVV equations |32|-|22|. 

2 Givental-style decomposition formulas for non-Gaussian CIV-DV 
partition functions 

• The shift of t-variables converts integral (|I.2.10|) into 



ZFw'it) = TFT" d^-^^ei^^ (H.2.1) 

Vc %(iV) J 

At the first glance, it can seem that there is no freedom in this definition and no place for the free 
parameters to get in. This is, however, incorrect, if one wants to allow any kind of analytical contin- 
uation in T-variables (coefficients of W(4>)). One of the ways to construct formal series (in powers 
of tk) solutions to the linear system of the Virasoro constraints (|I.3.10|) is provided by an analyti- 
cally continued saddle-point approximation. Different generators of the D-module are represented by 
expansions around different saddle-points <f> = (fro such that W'(4>q) = 0. If the polynomial 

n 

W'(x) = Y[(x- ai ) (II.2.2) 

8=1 
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has roots oti, then, since (fio are matrices defined modulo L r (A r )-conjugations (which allow one to 
diagonalize any matrix and permute its eigenvalues), the different saddle points are represented by 

4> = diag(ai, . . . ,«i;a 2 , . . . ,a 2 ; ••• ;a n , ...,a n ) (II. 2.3) 

with Oi appearing iVj times, Ys?=i Ni = N. 

Since we are dealing with a D-module, there is no need to keep these non-negative integers: in 
final expressions (like formulas for the multidensities and prepotentials) they can be substituted by 
any complex numbers. Moreover, iVj can depend on g and on (i.e. on the shape of W(<fi)), the only 
restriction coming from additional constraints (|I.3.16|) . The so restricted freedom to choose functions 
Ni(g,W) matches the free parameters / in expFg g.j. 

• The contribution of a particular extremum cj>Q (labeled by the set of Ni) is given by the Givental- 
style decomposition formula, 11 expressing it through the product of n Gaussian partition functions, 
each with its own Ni and M, = W"(ai) = ~~ a j) 



rr" c — NiW(a>i) -\7-_i n n n n 

4T 3 (#o) = n - ie Vnl „,„ Wo n «T Ni n °* n n z&^m ^ 

with operators 



Yol U(N) i<j " Kj i=l i=l 



6a = exp I 2 V (-) k (k+ k !r 1)] 4t 4t 1 > (H.2.6) 



k,i=o 



aj,- = aj — a,-, and 



*-(-s^) 

(W^( fc )(x) = 5*1^ (x)). The relation between the time-variables at the two sides of ea. (|11.2,5|) is given 
by 

oo/n \n/oo \ 

E t k E Tr *( a > + ^) fc = E E ( IL2 - 8 ) 

fc=0 Vi=l / i=l \fc=0 / 

with arbitrary JVj x JVj matrices <^>j. 

The operator V(z) = 2~Zfc>o IF+iW^ ls a nnear combination of the operators V(z)^ = 2~Zfc>o P^t - 1*7' 



V(^) =EV(^-«i)W (II.2.9) 



It is derived, as in by switching from an eigenvalue integral over N x N matrix (j> to n eigenvalue integrals 
over Ni x TVi matrices (f>i (each obtained with the shift by ac just changing of variables in integral 111.2. Il l, namely by 



rewriting the square of the Van- der- Monde determinant in (11.2.81 as 

N n I Ni \ n / Ni Nj \ 

= jj(^„ - = jj jj(^„ - ^, 6 ) 2 n n + - < ^) 2 = 

a<i> i=l \a<6 / i<j \a = l f>=l / 

= J] (&) n aJ iN * exp I 2 J] i ^b^ Tr ' Tr ^ ® ^ " J * ® ) = ( IL2 - 4 ) 

i — 1 \ m— 1 ^ / 

n n / oo , \ 

=n^(*>n«* 2 b-) 1 ^^ 

i=l i<j \ fc,i=0 »J ' ' / 

Despite the technical simplicity of this transformation, its real quantum-field-theory meaning (exact embedding of the 
SU(N) Cartanian r-function into the S (®™ =l U (Ni)) one) remains underinvestigated (in particular, no effective trans- 
formation in terms of matrices, avoiding intermediate eigenvalue integrals, is known), what makes generalizations to 
non-eigenvalue models obscure. 
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• Acting by operators V(z) on the decomposition formula (|II,2.5|) . one can express the non- 
Gaussian multidensities through the Gaussian ones. Note that V(z) = J2k>o z ~ k ^d/dt^ commutes 
with Oij and Oi and it does not act on Mj = W"{aj) = Y\j^i a iji which do not depend on t's. 
According to ()II.2.8|) and (|II.2,9|) . the action of V(z) on a particular Gaussian constituent, 

g 2 V{z)\ogZ^ U] {t^\N l )=pl h \z-a i ) (IL2.10) 

provides the corresponding Gaussian density, but with shifted arguments. As a result, we have, for 
example, for 

imatr) , , _ * M , yimatr) _ ™<J gg 6j (gU g (* = aj) gU 4^ (^g 

rii<j °ij Ui=l Oi [\i=l Z G ( tW \ N i) 

(II.2.11) 

Of course, Oi and Oij act both on Zq's and pc's. Similarly, expressions for the multidensities 
pw(%i, . . . , z m ) = V(zi) . . . V(z m ) log Zy/ contains as insertions in the numerator sums over split- 
tings of the set {zi, . . . , z m } into n subsets (perhaps, some of them empty, then the corresponding p 
with no arguments is just unity) 

E (f[p [ G i] (*n,-,^)) (IL2.12) 

{m}— >{mi),...,{m„} \i=l / 



• Expression (|II.2.11|) and its generalization with ()II.2.12j) are rather sophisticated. Role of the 
denominator is, as usual, to eliminate the disconnected Feynman-t'Hooft diagrams, including discon- 
nected contributions from different Gaussian constituents. The first contributions to (|II,2.11|) look 
as follows (we remind that pg( z ) = {{T(z\4>))) G , and introduce an obvious notation z% = z — on for 
shifted z-variables and 4>i for the N% x A^-matrix fields) 

p w (z) = p { ™ atr \z), 

P { w tr \z) = E ((rW^)e-^W))^ + E ((r^^Tl^e-^W))^ + . . 

3 (II.2.13) 

Here 

Wifc) = W(a t + fr) - W(a t ) - (H.2.14) 

and the structure of other corrections can be understood since every operator Oij inserts into the 
correlator a series 

00 / 2 \ 2 

exp (-) fc_1 7rX Tr i Tr i(^ " = 1 + —(NjTr^i - A^Tr^) + . . . (II.2.15) 

k=i \ ka ij J ai i 

The leading contributions to the l.h.s. and to the r.h.s. of the first line of (|II.2,13|) are respectively 



(see (11.3.441) and (11.3.461) ) and 



(W'(z))°- - iW'(z) YJ — ] =T.— + °(* 2 ) (II.2.16) 



E p [ £" M w - E f (* - f 1 - ^) - £ f + <W) (n.2.17) 
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(see (|II.1.13|) ). i.e. in the leading approximation (comp.with |3*9"| ) 



Si = gNi + 0{g 2 N 2 ) = Si + 0(S 2 ) (IL2.18) 

• The quantities p\y(z\t) and fw(z\t) are very similar to the non-trivial Gaussian averages, re- 
spectively, 

{{T(4>\z)))w = ((T(4>\z)e^ W ^)) G (II.2.19) 

and 

({Rw(<t>\z))) w = ((RwU>\z)*-** w )) (n.2.20) 

with 

1 W'(z) - W'(6) 

T^)=Tr -, R w (</>\z)=Tr U 7 W II.2.21 

z — <ft z — <p 

and W(z) = W(z) — W(a) — ^-(z — a) 2 . The problem is, however, that these averages are not 
well-defined, because there are different saddle points a, which are roots of W'(z) (since these are 
saddle points, they can be both minima and maxima of W(z)), with different values of M = W"(a). 
Expressions (|11.2.19|) - (|11.2.20|) can make sense only for particular W, associated with one of the roots. 
Thus, they describe only particular, perturbative branches of the partition function, when just a single 
Ni = N and all other Nj = 0. To describe generic non-perturbative branches, one needs to return to 
(111.2.111) and (111.2.121) . 

As to perturbative phases, in any of them one can take the expansion pi. 1.3(1 for the Gaussian 
prepotential, 

J$ = F» + j vp^ + - j j vvp^ + ... (II.2.22) 

and define an expansion for the non-Gaussian one, by simply performing a shift v(z) — > —W(z) +v(z), 
i.e. by resummation of the diagrams 

so that 



J%> = F^+j> vp^ l) + -J <f> vvp^> + ... (II.2.23) 



Fw = F G + jwp^ + <fwWp { f ] + ... , 
# ) =4 |1) +/* + y/^P? + - : 

OO -. OO -. 

k=0 k=0 
J fci,fc 2 =0 Z l Z 2 



(II.2.24) 



3 CIV-DV prepotential 



As we saw in part I, for a given W the solutions to Virasoro constraints, which possess the genus 
expansion, are labeled by arbitrary function F[g,T] of the first n time variables To, . . . ,T n -± and g 
(dependencies on T n and T/v+i are prescribed by p.3.43Jl ). To match the generic F[g, T], the variables 
Si in the decomposition formulas pi. 2. 5(1 should also be made T- and g-dependent (though their 
possible dependencies on T and g are severely constrained). 
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The CIV-DV prepotential |13]-[37] is defined as F[g,T] with constant 



Si = const (II. 3.1) 

Then it is immediately given by (111.2.51) . but it remains unclear, what kind of additional constraint 
should be added to the Virasoro equations 1)1.3. 10|) in order to reproduce this condition 1)11.3.1)1 . Instead, 
as was discovered by R.Dijkgraaf and C.Vafa a new representation occurs, that is, in terms 
of Seiberg-Witten theory PS|,|3Hl-[52- See Tables below for explicit expressions for the CIV-DV 
prepotential 

F DV [T\S] = Y J 9 2p ^ ) y[T\S] (H.3.2) 

p>0 

• In the framework of the Seiberg-Witten-Dijkgraaf-Vafa theory the T-independent S% can be 
expressed through J^y 



Si = f yw(z)dz 
J A, 

yl w ( Z )=[W'(z)] 2 -4f° V (z), 

r(°) fT-ICl _ ST fl. , ; , o\rr „k OJ ~t>y 



f^{z) = R w {z)T^ v [T\S] = £ (k + l + 2)T k+w z k9:F >> : 



OF, 



(0) 



DV 



dS„ 



Bi 



yw{z)dz 



Now we are going to calculate the CIV-DV prepotentials using the decomposition formulas 1)11.2.5)1 . 
• We begin with defining the action of the operator 

a, = ex P (- £ = exp (- M ;t (.1.3.4) 



with 



n 



Mi = W^(oi) = W"{a t ) = {n + l)T n+l J[ a 
2M l a^ ) = W^(a l ) = 2M l J2- 



6M^ = W^( ai ) = 6M t J2 (H.3.5) 



(p + l)\M i( jf = W^ +2 \ ai ) = (p + l)\Mi J2 



1 



31,— Jp¥=i 
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on a particular Gaussian partition function ZQ l {t\Ni): 

( 

= 1 + 

t=o 



O l Z^{t\N l 



V 



2!(3!) 2 



(W 3 ) 



+ 



W^(«0 



6! 



M, 



3!5! 



tr^ 3 tr<^ 



Mi 



+ 



^(4) ( 



a. 



G 2!(4!) 2 



W^( ai ) W^( ai ) I z 
2W4! ^ 



1 



+ 



G 4!(3!)4 



2 \Mi 
I G 

(tr^ 4 ) 
\ 



+ 

2\ M i 

G 



Mi 



G 



+ 



+ 



W( 3 \ai)W( 7 \a>i 



8! 



t^ 8 ) G + 



3!7! 



trc/r tT(/>' ) + 



W^iaAW^iaA I j4 , fi \M 



^(5)(. 



2!(5!) 2 



(tr0 5 ) 



2 \Mi 
G 



2!(3!) 2 4! 



(tr(p 3 y tr0 £ 



M, 



3!(4! 



H3 



(tr^ 4 ) ; 



+ 



3!(3!) 3 5! 



3 \ M 4 

5 



(tr0 3 ) tr<£' 



+ 



G 



z%*{t\Ni) 



t=0 



(II.3.6) 



• This formal series becomes ill-defined when some ctij = 0, i.e. when we consider an expansion 
near the degenerate extremum aij = ctj. In such cases (for such W(z)) one should develop a different, 
essentially non-Gaussian formalism, involving expressions like (for N = l) 12 



/ #ex P (-w + £ t k A = y: ^# 

C \ k=l I k>0 p h p 



C{ P ) 



k + 1 

P 



(II.3.7) 



where Shk{t) are Shur polynomials, 



exp =E S/l *W* 

\fc=l / fc=0 



and the Gamma-function for the chain C is defined as 



Y c {z) = I (ff'^e-*^ 
Jc 



(II.3.8) 



(II.3.9) 



Generalization of this calculation for N ^ 1 is essentially provided by the non-Gaussian multiden- 
sities p for W(z) = W p (z) = hz p . Note only that association of chains C with the function F\y and 
the polynomials fw p is somewhat non-trivial: for example, a naive choice of /C 3 ! 1 ) = u = const would 
imply that the correlators 

tr^) (0) = I z k p§ l \z\f = const)dz = - \ I z k Jp 2 z 2 P~ 2 - ASdz (II.3.10) 

' J CO ' £ J oo 

do not vanish only if k + 1 is an odd multiple of p — 1, i.e. if k = (2n — l)(p — 1) — 1, while for 
N = 1 it should be equal to Tc (^r)- (I n the Gaussian case p = 2, this requirement for C is that 
12 The contours in the sum depend on p due to the change of integration variable which depends on p. 
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Tq(z) 7^ only for half-integer, but not for integer arguments.) Note also that, for p 7^ 2, the choice 
of /(oil) 

= v = const is in no way distinguished, since there is no M and thus no natural definition for 
v (which is equal to gN/M in the Gaussian case). 

In what follows we consider the non-degenerate situation, when all ctij 7^ 0. 



Even if one needs the decomposition formula (|II.2.5f) at t = 0, the terms like 



3! dt 3 



m 



6w(cti), which do not contribute to (jll.3.6f) . should be kept, since they can combine with Oij and give 
non- vanishing contributions even at t = 0. 

For explicit calculations, we introduce a double gradation of the prepotentials: by powers of g and 

S, 

.F = E F {P) = E ^ (P ' m) + T = E Frn + F 



p>0 



p>0 



m>l 



(II.3.11) 



where J 7 (pl m ) and T m denote contributions with the m-th powers of S"s of genus p and of all genera re- 
spectively. T is the contribution singular in S, which can be extracted from II. 2^91 Explicit expressions 
(with finitely many terms) can be found for every particular J r ^'\ m ) . The decomposition formulas 

pi.2.5|) express T\ 



(p\m) 
DV 



through Tq ' m ^ with p' < p, m' < m, and some t 7^ 0: at the end, one 



t=o 



t=0 



should put the lowest times 



^4-2 

k • 



3 < k < n+ 1 



(II.3.12) 



but, before doing this, one should take some ^derivatives. 
The decomposition formula ()II.2.5|) states that 



-(Sd)n 2 (cdd)(l/g 2 )F G 



where the "normalization" part of the prepotential 



e (i/g 2 )F p D 7 1 



n 



II, Vol U{Si/g) Y M ^/2 9 : 
the sum of the Gaussian prepotentials 

^=e^=x>g i (* w ) 

m>l i=l 

and the operators 



n 



2 11 



(II.3.13) 



(II.3.14) 



(II.3.15) 



ka k - a+W 



(cdd) = e E 

i+i \k,l>l 



(II.3.16) 



(-) k+1 (k + 1-1)1 d d_ 



Note that (c<99) annihilates any single -F m , it acts non-trivially only on bilinear and higher-order 
combinations of F's. 

Substituting (see [2H] ) 



log (Vol U{S/g)) = log(|) + 1 log(|) + S " I; l0 S( 2 ^ " C'(-l) " E ^ ^ 



(J 



1 

12 



, 5* N 
9 



3ff 2 5 
V ~ 2^ 



2fc-2 



.^ 2 4fc(fc - 1) S 2k ~ 2 
(II.3.17) 
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one straightforwardly obtains from (jll.3. 13|) (by expanding exponentials at both sides and comparing 
the terms with the same powers of S and g 2 ) 



DV 



t=Q 



DV 



\t=0 



J2SiW( ai ) + Fi| t 



k k 



■i - 4S ^ + 5 D lo s - 1 E 



1 00 2 fe-2 

F 2 -(5a) J p 1 + -^^ r ( c 9a) fc F 1 2 
fe>i 



k k 



3 



t=0 



i 71 3-(^)F 2 + -(5a) 2 F 1 + 



00 „2fc-4 



00 2k-2 1 * : 

■ E ^Tr(cdd) k [FiF 2 - F!(Sd)Fi] + - E ^7T- ( - ZF^cddfF 2 
fc>i k[ 6 fc>2 fc! V 



k k 



(II.3.18) 
, with 



More terms of the expansion can be found in the Tables. 

In the sums still remaining in (111.3.181) . only finitely many terms contribute to J-^y 

any given p and m, and one should also take into account the non-trivial g 2 -dependence of F n 
ZiF m (Si\t®) 



F m {S\t) = S m Y J 9 2p 

p>0 



OO OO jil jls 

k\ ks 



E E 

S = l k 1 <...<k s 
ll,...,l s 



h\...l s \ 



k l1 h ls 

! • • • ) 



(p) 



(n.3.19) 



M2 



-j S\- i yUki-2)=m + 2p-2 

±(hk l +...i s k s ) \ 2 4^ x n ' 



where we introduced a condensed notation for combinatorial quantities (calculating the numbers of 
the relevant fat graphs) 



(k[\..., k l f)) iP) = M ^+...l s k s ) //(tr^ 1 . . . (tr«^) 



G 



(II.3.20) 



The delta-function in 1)11.3. 19|l restricts the possible number of non-vanishing = (with 3 < k < 
n + 1) for given m,p, t\ and i 2 . However, since at the end of the day, we are going to put t\ = and 
t<i = 0, they should be eliminated by applying the operators (Sd) and (cdd), which contribute extra 
powers of 5 and g 2 respectively. Thus, only finitely many terms in the expansion (111.3.19)1 contribute 
to a given P^v . A useful way to study these selection rules is to introduce the notation for 
lj in 1)11.3.2011 . Then, the delta-function in (111.3. 1911 imposes the constraint 



for the correlator 



to be non- vanishing. 



E(fc - 2 Wk = M + 2m + 4p - 4 

k=3 

(({M,AA 2 ,...})) = (((tr^(tr0 2 ) A '' 2 ...)) 



(II.3.21) 



(II.3.22) 
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• Particular example of m = 1 : 
Since 



DV 



+ J2s i W(a i ) = Ft\ ii)= _ T(i) 



(II.3.23) 
(II.3.24) 



and in (111.3. 19j) 

CO 

^(fc-2)A4=M + 4p-2 

fc=3 

there are no non- vanishing contributions to the r.h.s. for p = 0, while for p = 1 we have 2A/4 + A3 = 2 
(i.e. either A/4 = 1 or A/3 = 2), and for p = 2 we have 6A/g + 5A/V + 4A/" 6 + 3A/" 5 + 2A/4 + A/3 = 6: 



(o|i) 



DV 



t=0 



£S^( 



a; 



T 
•1 i 



DV 



S, 



t=Q 



y*. 

^ M 



U) , ((3 2 )) (1) 
4 ^ + 2!-3 2 



(7 



(II.3.25) 



5; 



y*. 

4- Mi 



C7 (i) 
°2 



+ 



Si 



2! -3 2 



12M- 



,0 



3(7. 



(0 

2 J ' 



See Tables for extra explicit expressions. 
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Part III 

Tables 



1 Gaussian partition function 

1.1 First generalized Catalan numbers 

The number of diagrams of different genera in the "leaf graph (Tr<j) 2k ) 

2k total# p 1 2 3 4 5 

(2k-l)\\ 






1 


1 












2 


1 


1 












4 


3 


2 


1 










6 


15 


5 


10 










8 


105 


14 


70 


21 








10 


945 


42 


420 


483 








12 


945 • 11 


132 


2310 


6468 


1485 






14 


945 -11-13 


429 


12012 


66066 


56628 






16 


945 -11-13-15 


1430 


60060 


570570 


1169740 


225225 




18 


945 • 11 • 13 • 15 • 17 


4862 


291720 


4390386 


17454580 


12317877 




20 


945 • 11 • 13 • 15 • 17- 19 


16796 


1385670 


31039008 


211083730 


351683046 


59520825 



(III.1.1) 



1.2 First Gaussian densities (generating functions of the generalized Catalan num- 
bers); v = gN/M 



Genus p = 0: 



where 



• Genus p = 1: 



M P [) ~ 2 ~ ^ 



4 fc r(A: + l/2) u k+1 

^ r(i/2)(fc + i)!^+i 



y 2 {z) = y 2 G {z) = z 2 -Au, 
gN S 



v = 



M M 
^ A k T{k + 5/2) u 2k+1 



yHz) r ( 5 /2)£! 



(III.1.2) 



(III.1.3) 



(III.1.4) 
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• Genus p = 2: 

3 (2|1) 21«/ (z 2 + »j ~ / 4 fe r(fc + ll/2) 4^1^ + 9/2) \ 

^ 1 J „iu^ Z1 Z^ rm/oM-i r rm /out _ im ^2 



L W fbA r(n/2)jfe! r(n/2)(jfe-i)! / z 2fc + 9 

° (III.1.5) 

y 14^4 fe - 1 r(A; + 9/2)(5fe + 18)i/ 2fe + 1 

z» + Y ^ r(9/2)(fc-l)! ^fcT9 



• Genus p = 3: 



M 5 (3|i) (;s) = 11^(135z 4 + 558^ 2 + 158^ 2 ) = 

y 17 (z) 

^/ 4 fe r(fc + 17/2) 4^1^ + 15/2) 4 fc - 2 r(fc + 13/2) \ 

11 2. I W5 ■ r(17/2)fc! + 558 ' r(17/2)(ife - 1)! + 58 ' r(17/2)(ife - 2)! J ^+13 



fc=0 



v v 



2 



+ v 15 + 1 Q . 1 K X] 



44 ~ 4 fc - 2 r(fe + 13/2) (aA; 2 + 6A; + c) z^ 2fc+1 



z 13 z 15 13-15^ r(13/2)(fc- 1)! z 2fe + 13 

(III.1.6) 

ak 2 + bk + c= 135 • 16(fc + 13/2)(fc + 15/2) + 558 • 4k(k + 15/2) + 158fc(A; - 1) (III.1.7) 



• Generic genus p > 0: 

In general, 



(III.1.8) 



M 2p-i . = I") ( „ LL9) 



where 



(III.1.10) 



is a polynomial of degree p — 1 in z 2 , 

Q Q (x 2 \u) = v 
Qi(x 2 \u) =21(x 2 + u) 
Q 2 (x 2 \v) = ll(135a; 4 + 558x 2 z/ + 158z^ 2 ) 
Q 3 (x 2 \u) = 143(1575x 6 + 13689x 4 z^ + 18378xV + 2339z/ 3 ) (III.l.ll) 
Q A {x 2 \v) = 88179(675x 8 + 9660x 6 z/ + 28764xV + 18908xV + 1354z/ 4 ) 
Q 5 (x 2 \u) = 111435(218295x 10 + 4534875x 8 z/ + 23001156xV + 34604118xV 

+13447818xV + 617926z/ 5 ) 

1.3 First recurrent relations for Gaussian mu/tzdensities 

Some recurrent relations for the Gaussian multidensities 

y W^fa.«,)- ^" w -' W " w (in.i.12) 

Z\ — Z2 

y(z 1 )pW 1 \z 1 ) = pM 2 \z 1 ,z 1 ) (III.1.13) 

(III.1.14) 
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y (z 1 )p^\z 1 ,z 2 ,z 3 ) = 2p(°\ 2 \z 1 ,z 2 )p^ 2 \z 1 ,z 3 ) + 

p^ 2 \ Zl ,z 3 )-pm(z 2 ,z 3 ) , a p("l 2 )(z 1 ,z 2 )-/}(°l 2 )(z 2 ,z 3 ) 
z\ - z 2 Zi- z 3 



(III.1.15) 



y(z 1 )p^Hz 1 ,z 2 )=2p^Hz 1 ,z 2 )p^\z 1 )+^ (IIL1.16) 

z\ — z 2 

y( Zl )p^( Zl ) = (p^\z 1 )) 2 + p^(z 1 ,z 1 ) (IIL1.17) 

(III.1.18) 

1.4 First Gaussian mu/ttdensities: analytic expressions; z/ = gN/M; y(z) = z 2 — Av 

The first Gaussian miittidensities (one-, two-, three-, . . .point functions) are: 
• Genus p = 0: 

_L p <oii) (2) = iz|(£i (m.i.19) 

^-^^(OT- 1 ) (ra - L20) 

there is no singularity at zi = 22, for coinciding variables the genus zero two-point function 

pm ^ Z) = ^{z) ( IIL1 - 21 ) 

M/> (0|3) (zi z 2 z 3 ) = — ( - + 

y(zi)y(z 2 )y(z 3 ) \( Zl - z 2 )( Zl - z 3 )y 2 ( Zl ) (m 1 22) 

1 1 \ _ 1v(z\z 2 + z 2 z 3 + z 3 z\ + Au) 



+ 



- ^i)(^2 - z 3 )y 2 {z 2 ) (z 3 - z 1 )(z 3 - z 2 )y 2 (z 3 )J y 3 (zi)y 3 (z 2 )y 3 (z 3 ) 

*' V "" ( '* M) = R^fej^SW * 

X {^\Az\z\z\z\{z\Z 2 + Z1Z3 + Z1Z4 + Z 2 ^3 + -^2-24 + Z 3 Z±) 

+3ziz 2 z 3 Z4(zfz 2 z 3 + z\z\z\ + z[z\z\ + z 2 z 3 z\)\ 

-8[(zfz 2 Z 3 + z\z\z\ + z\z\z\ + Z 2 Z 3 z\)(ziZ 2 + + Z1Z4 + + Z 2 2 4 + Z3Z4) 

1 c ^ 2 2, 22, 22, 22, 22, 2 2\ c 3 3 3 3i 

+QZ\Z2ZzZ\{z l Z 2 + + z l z 4 + + z 2 z A + 2 3 z 4 ) - 6z 1 z 2 ^3^4jz / 

-32[2(zfz 2 + Z 2 z| + Z 2 Z 2 + z\z\ + Z 2 ^4 + z\z\){z\Z 2 + Z1Z3 + + Z 2 Z 3 + Z 2 ZA + z 3 za) 

+3(zfz 2 z 2 + z\z\z\ + z\z\z\ + z 2 z 3 z\) - §z\z 2 z 3 Zt±(z\ + z 2 + z 2 + z\)\v 2 

+ 128[5(z 2 + z\ + Z 3 + Z4)(ziZ2 + Z1Z3 + ^1^4 + -^2-^3 + -^2-^4 + Z 3 Za) - \2z\Z 2 Z 3 Za\v 3 

+512[3(z 2 + z\ + z\ + z|) - 8(21 z 2 + z x z 3 + 212:4 + z 2 z 3 + z 2 z± + z 3 Zi)]v A - 12288i/ 5 ' 

(III.1.23) 

• Genus p = 1: 

Mp(1|1)( " ) = ?M ( m - 1 - 24 ) 

M 2 ^ 1 ! 2 ) (2:1,2:2) = 77 ^77 , ^i^(5^ + 4zfz 2 + 3z 2 z| + 4ziz| + 54)+ 

y (zi)y (z 2 ) V - (HI.1.25) 

+4^ zf - 13ziz 2 (z 2 + ziz 2 + zf ) + 4 + 16^ 2 (-z 2 + 13ziz 2 - z|) + 320^ 3 ) 
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" y{*ify{*2?y{z 3 f 

{^z\z\zl{z\z\ + z\zl + z\z\ + z x z 2 z 3 {z\z\ + z\4 + 44 + *i*s + z\z\ + 44)) 
+\z\z\z\(z\z\ + zfz 2 + zizf + zfz 3 + z 2 ^f + z|z 3 + z\z\ + 44 + z\z\) 
+2>z\zlzl{ Zl +z 2 + z 3 )] + i[5(z lZ2 + ZiZ 3 + Z2Z 3 )(44 + 44 + 44) 

-ios444(44 + 44 + 44 + 44 + 44 + 44+ 

Z lZ \ Z \ + Zi44 + 4^4 + z\z 2 z\ + 44*3 + 4z2 z 3) 

-96zfz% z\{z x z 2 zl + zizfz 3 + z\z 2 z 3 + + z\z\ + zfzg + z?z§ + 44 + 44) 
-\8foz\z\z\(z\Z 2 + Z\Z 3 + Z 2 Z 3 ) + Yl44*\{*\*\ + 44 + *2*i)\ v 

+ie[ - (44 + 44 + 44) - 21(44(4 + 4) + 44(4 + 4) + 44(4 + 4) 

+z 1 z 2 z z (44 + z\z\ + 44 + z\z\ + 44 + 44)) 

-27zi z 2 ^(44 + 44 + 44 + 44 + 44 + 44) 
~^444(44 + z ?4 + 44 + z ?4 + 44 + 44) + 1324222! (4 + 4 + 4) 

+\^z\z\z\{zxz\ + z\z 2 + 2i zf + zfz 3 + z 2 4 + 4*3) 
+3784 z\z\{z\z\ + z\z 2 + z\z\ + 423 + z 2 z\ + 423) 
+?&<oz\z\zl(z\zl + afzjj + zfzf) + 234^44]^ 

+i92[ - (44(4 + 4) + 44(4 + 4) + 44(4 + 4)) 
+60444(44 + 44 + 44) + 9(44(4 + 4) + 44(4 + 4) + 44(4 + 4) 

+z\Z 2 z 3 (z\z\ + z\z 2 + z\z\ + z\z 3 + z 2 z\ + z\z 3 )) 

+n(z\z\ + 44 + 44) - 156444(^4 + 4*2 + *i*l + 4^3 + z 2 4 + 4*s) 
+48444(4 + 4 + 4) + i3ziz 2 z 3 (44 + 44 + 44 + 44 + 44 + 44) 

-144444(^1 + 22 + -23) + 3iziZ2z 3 (44 + 44 + 44 + 44 + 44 + 44)]^ 3 
-256[3(44 + 44 + 44) + 6(44(4 + 4) + 44(4 + 4) + 44(4 + 4)) 

+129(44(4 + 4) + 44(4 + 4) + 44(4 + 4)) 
+52z!z 2 z 3 (4 + 4 + 4) + 46(ziz 2 (4 + 4) + 2123(4 + 4) + W4 + 4)) 

+378444(21 2 2 + ziz 3 + z 2 z 3 ) + Yb?>zxz 2 z 3 (z\z\ + z\z\ + z\z\ + z\z\ + z\z\ + z\z\) 
+ 'bAQz\z\z\( K z\ + 4 + 4) + y^z\z 2 z 3 (z\z 2 + 4z 2 + *1*3 + 4*3 + *2* 3 + z 2 z 3 )]v A 

+io24[ - 10(4 + 4 + 4) + 75(44(4 + 4) + 44(4 + 4) + 44(4 + 4)) 
+i74( Zl z 2 (4 + 4) + ,^(4 + 4) + z 2 2 3 (4 + 4)) + 237(44 + 44 + 44) 

249z 1 z 2 z 3 (zfz 2 + zi4 + 4z 3 + *i*3 + 4*3 + 224) + 186ziz 2 z 3 (4 + 4 + 2!) + 828zfz 2 2 z 3 ]v l 
-12288[80(ziz 2 (4 + 2!) + 2iz 3 (4 + 4) + 2 2 2 3 (4 + 4)) - 12(2? + 4 + 2!)+ 

(III. 
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+39(z 2 z 2 , + z\z\ + z\zf) + 82212:223(^1 +z 2 + z 3 )]v 6 

(III.1.27) 

+344064[ - 2(2? + z\ + z\) + 7(ziz 2 + ziz 3 + 2; 2 -Z3)]^ 7 + 3342336z/ % 



• Genus p = 2: 



M 3 p (2|l) (z) = 2M^ + ^ (mi28) 



M 4 /9 (2|2) (zi,z 2 ) = — — ^— —-(i9Qz^4 + 16&4°4 + 147z\z\ + ^44 + 165x14+ 
Hi--:)' '!Ji '■■>'■■ V - - 

i5Gzf4 + 14744 + 1^44° + i^44 l 

+ [399^ 1 2!2 + 462z\°z% ~ M1144 ~ 3664zf z\ - 331744 ~ 334244 
-33\7z\z\ - 3664zfzf - 4011^4 + 462zf 4° + 399zi4V 
+big[84zl° - 9408z?z 2 - W516zfz^ + 33736zja| + 30980^4 + 27624zf z\ 
+30980z 4 4 + 33736z?4 - 10516z?zf - 9408zizf + %4zf\v 2 
+ [ - \m<oz\ + 92368z[z 2 + 101120*?jzf - 129648^1 - 133248z 4 z 4 
-129648z 3 zf + \§\\2§z\z\ + 92368zi4 - 1696zf]i/ 3 
-[12992z? - 487936*? z 2 - 491072z 4 z\ + 202752^ z\ - 491072z 2 4 - 487936*i*| + \2m2z\\v^ 
+ [ - 7731224 + 1415936z 3 z 2 + 2519552z 2 z| + 1415936^1 z\ - 77312zf]i/ 5 
+ [ - 560128z 2 - 2951168ziz 2 - 560128z|]i/ 6 - 1736704i/ N 



(III.1.29) 

• Genus p = 3: 

M5 (3]1) = 11,(135^ + 558^ + 158^) (HI. 1.30) 

Note that all these multidensities are defined for any finite value of N and have nothing to do with 
the limit N — > 00. See also |53| for an alternative approach to the similar problem. 

1.5 First Gaussian mu/tadensities: first terms of iV-expansions 

In this section M = 1 and g = 1 so that v = N. 

• One-point functions: 

p®%) = ±N + liV 2 + V 3 + V + + ^ + ^ + ^ + ^ 9 + °( Nl °) 

rf* rf* O /y* O /y* I nf \j rf* J- _1_ rf* -L O rf* -L O rjf* J- I 

(III.1.31) 



t*J tti iA/ <Jy •*> 



(HI)/ X 1 »r 10 *r2 70 ,.3 420 , 2310 . 

pg |1J x = -iV + -iV 2 + -iV 3 + — iV 4 + — JV 5 

lL it/ lL lL 

12012, 60060 at7 291720 atS 1385670 q ^ /Arlt ^ 
+— N + — Nl + + + 0( ~ N ) 

(211), s 21 r 483 at2 6468 Ar , 66066 Ar4 570570 Ar . 
4390386 ATfi 31039008 at7 205633428 at8 1293938646 q 

AT® i AT' I AT° l AT-- 1 



j.^rrp:^ + """"""" ^ + — tv« + + o(iv 10 ) 



/yjl_9 /yj21 /yj23 /yj25 

(311),, 1485 Ar 56628 Ar2 1169740 , 17454580 Ar4 211083730 

Pa ( x ) = -zvrN + ~^r^ + — i^iv 3 + ^r^^ 4 + zii 

*X' iL iL 



iV 



5 



2198596400 ATfi 20465052608 at7 174437377400 , 1384928666550 atQ ^ /ATin , 
== A^ 6 + = N + = N 8 + iV 9 + 0(N 10 ) 

^Jj *)C *jC 
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(III.1.32) 



(III.1.33) 



(III.1.34) 



• Two-point functions: 



f |2), . 1 f AN\ ( 2N QN 2 20iV 3 70iV 4 252iV 5 \ 



N 

X 2y2 



+ N 2 



+N 



+N 5 



126 



+ 



+ 




X 9y3 



+ 0(iV 6 ) 

(III.1.35) 



■*x v y 

70 60 60 



x Sy2 



x 7y3 



+ 



630 

x 10y2 



560 630 600 600 



x 9 y 3 



+ 



x 8y4 



X 7y5 



x 6y6 



600 630 
+ 



x 5 y 7 



x 4 y 8 



+ 



560 

X 3y9 



+ 



630 

X 2yl0 



)n 3 



+ 



4620 4200 5040 4760 4900 4800 



+ 



+ 



x 10 y 4 



+ 



X 9y5 



+ 



X Sy6 



x 7 y 7 



+ 



4900 4760 5040 4200 4620 



+ 



+ 



+ 



+ 



x 2 y 12 



N 



/ 30030 

V x 14: y 2 



+ 



x 6y8 x 5 y 9 x 4 y 10 x 3 ?/ 11 
27720 34650 32760 34650 33600 34300 



x 13 y 3 



+ 



12 4 

x y 



+ 



+ 



x io y 6 



+ 



x 9 y 7 



+ 



y°x° 



+ 



+■ 



33600 34650 32760 34650 27720 30030 



x 7 y 9 



+ 



X 6yl0 



X 5yH 



pg' 2 W) = (^ + X 8,2 



189 168 



+ 



147 

x 7 y 3 



+ 



+ 



x 4 y 12 



+ 



x 3 y 13 



+ 



x 2 y 14 



N 5 + 0(7V 6 



156 



+ 



165 

x 5y5 



156 

x 4 y 6 



+ 



147 

x 3 y 7 



+ 



168 189 



X 2y8 



+ —) N 

x"y 



+ ( 



5313 4830 4725 4760 4795 4770 



x u y 



+ 



x 10 y 2 



+ 



X 9y3 



x 8 y 4 



+ 



x 7y5 



x 6y6 



4795 4760 4725 



x 5 y 7 



+ 



X 4y8 



X 3y9 



4830 5313 , Ar2 



X 2yW 



xy 



11 



+ ( 



84084 77616 81774 80640 80640 80640 



X l3y 



x 12 y 2 



+ 



+ 



x 10 y 4 



+ 



X 9y5 



+ 



X Sy6 



81340 80640 80640 80640 81774 77616 84084 



re 7 ?/ 7 



+ 



X 6y8 



X 5y9 



+ 



+ 



X 3yH 



x 2 y 12 



+ 



xy 



13 



)iV 3 + 0(iV 4 ) 



(III.1.36) 



(III.1.37) 



Three-point functions: 



p^ 3 \x,y,z) = ( 



2 2 2 . 

+ n 5 o + n „ „ 1-ZV 



x 2 y 2 z 3 x 2 y 3 z 2 x 3 y 2 z 2 



+ 



12 



+ 



12 



+ 



12 



x 2 y 2 z 5 x 2 y 5 z 2 x 5 y 2 z 2 



+ 



12 



+ 



12 



+ 



12 



+ 



12 



+ 



12 



+ 



x 2 y 3 z 4 x 3 y 4 z 2 x 4 y 3 z 2 
12 



+ 



+ 



x 2 y 4 z 3 x 3 y 2 z 4 x 4 y 2 z 3 x 3 y 3 z 3 < 
60 60 60 : 60 

p3/i/6-y2 /-yi6/ij3-y2 <y6n"3/j|2 -y3^y>2/i<6 -y3 /y»6/j i2 <y6/y>2/i|3 



60 



+ 



60 



+ 



-x*'y y 'z^ x"y ,J z* z^x^y" z^'x^y" z^x^y* z"x^y u 



(III.1.38) 



36 



60 

x 7 y 2 z 2 
72 
x 5 y 4 z 2 



+ 



+ 



60 

z 7 x 2 y 2 
72 



+ 



60 

x 2 y 7 z 2 
72 



+ 



z 5 x 4 y 2 



+ 



Z i X 5y2 



72 

z 3 x 4 y 4 
72 



+ 



72 

z 4 x 4 y 3 
72 



+ 



z 4 x 2 y b 



+ 



72 

z 4 x 3 y 4 
72 



48 



z 3 x 3y5 



+ 



48 



48 



z 5 x 3y3 



z 3 x 5y3 



+ 



280 



+ 



280 



+ 



280 



+ 



280 



+ 



280 



+ 



z 5 x 2 y 4 x 4 y 5 z 2 
280 280 280 



+ 



+ 



+ 



280 



z 9 x 2y2 x 9y2 z 2 x 2y% z 2 z 3 x 8 y 2 z 8 x 3 y 2 z 3 x 2 y 8 z 8 x 2 y 3 x 8 y 3 z 2 x 3 y 8 z 2 
360 360 360 360 360 360 



z 4 x 7 y 2 



360 



z 5 x 6y2 



360 



z 4 x 6 y 3 



+ 



+ 



+ 



z 7 x 4 y 2 
360 

z 6 x 5y2 

360 

z 3 x 6 y 4 



+ 



+ 



+ 



z 4 x 2 y 7 
360 

z 5 x 2y6 

360 

z 6 x 4 y 3 



+ 



+ 



+ 



z 7 x 2 y 4 
360 

z 6 x 2y5 

360 

z 3 x 4 y 6 



+ 



+ 



+ 



x 7 y 4 z 2 
360 

x 5y6 z 2 

360 

z 6 x 3 y 4 



+ 



+ 



+ 



x 4 y 7 z 2 
360 

x 6y5 z 2 

360 

z 4 x 3 y 6 



+ 



240 

z 3 x 7y3 



+ 



240 



+ 



240 

z 7 x 3 y 3 



+ 



288 

z 3 x byb 



+ 



288 



+ 



+ 



1680 

\z 3 x 8 y 4 

1680 
x 4 y 9 z 2 
1680 

z 8 x 5y2 

1440 

z 3 x 7 y 5 
2160 

z 4 x 5 y 6 



+ 



+ 



+ 



2160 

' z 4 x 7 y 4 
1800 

z 7 x 6y2 



+ 



1680 

z 4 x 8 y 3 
1680 
z 9 x 4 y 2 
1680 

z 5 x 2y8 

1440 

z 5 x 7y3 

2160 

z 6 x 5 y 4 
2160 
z 4 x 4 y 7 
1800 

z 6 x 7y2 



+ 



+ 



+ 



1680 

z 4 x 3 y i 
1680 
z 4 x 2 y 9 
1680 

z 5 x 8y2 

1440 

z 5 x 3 y 7 
2160 

z 5 x 6y4 



2160 

z 7 x 4 y 4 



+ 



+ 



+ 



+ 



+ 



+ 



+ 



+ 



+ 



+ 



288 

z 5 x 5y3 



1680 

z 8 x 3 y 4 
1680 
z 9 x 2 y 4 
1680 

z 8 x 2y5 

1440 

z 7 x 3 y 5 
2160 
z 6 x 4 y 5 
1120 

z 3 x 3 y 9 



+ 



+ 



+ 



+ 



+ 



+ 



+ 



432 

z 4 x 5 y 4 



1680 

z 8 x 4 y 3 
1680 
x 9 y 4 z 2 
1680 
x 8 y 5 z 2 
1440 

z 3 x 5y7 

2160 

z 4 x 6 y Fj 
1120 

z 9 x 3y3 



+ 



432 

z 4 x 4 y 5 



+ 



432 

z 5 x 4 y 4 



)n 4 



+ 



+ 



+ 



+ 



+ 



+ 



+ 



1260 



z 3 x V)y2 



+ 



1800 

z 6 x 3y6 

1260 

z l0 x 3 y 2 



+ 



+ 



1800 

z 3 x 6y6 

1260 

z 3 x 2 y 10 



+ 



+ 



1800 

x 7y6 z 2 

1800 

z 6 x 6 y 3 
1260 
z 10 x 2 y 3 



+ 



+ 



+ 



1800 

z 7 x 2 y 6 
1260 
z n x 2 y 2 
1260 

x 10 y 3 z 2 



+ 



+ 



+ 



1800 

z 6 x 2y7 

1260 

x n y 2 z 2 
1260 
x 3 y 10 z 2 



+ 



1680 

z 3 x 4 y 8 
1680 
z 4 x 9 y 2 
1680 

x 5y8 z 2 

1440 

z 7 x 5y3 

2160 

z 5 x 4 y 6 
1120 
z 3 x 9 y 3 
1800 

x 6 y 7 z 2 



+ 



+ 



+ 



1260 

x 2 y n z 2 
1728 

z 5 x 5y5 



P ( G 3 \x,y,z) = ( 



24 



x 3 zy 4 



+ 



24 



24 



z 2 x 2 y 4 



+ 



24 



z 2 x 4 y 2 



+ 



30 



x 2 zy 5 



+ 



30 



z 5 x 2 y 



+ 



30 



z 2 xy 5 



+ 



x 4 zy 3 

24 
z 4 x 2 y 2 
30 

+ ■ 



+ 



24 



z 4 xy 3 



+ 



24 



z 3 x 4 y 



+ 



24 



z 4 x 3 y 



+ 



24 



N 5 + 0{N 6 ) 

(III.1.39) 



z 3 xy 4 



+ 
30 



18 



z 2 x 3 y 3 



+ 



18 



z 3 x 2 y 3 



z 5 xy 2 



Z 2 X 5 1 



+ 



30 

x 5 zy 2 



+ 



30 

x 6 zy 



+ 



+ 



18 



z 3 x 3 y 2 



30 



+ 



+( 



540 



x 5 zy 4 



+ 



540 

z 4 x 5 y 



+ 



480 



z 2 x 5 y 3 



+ 



540 



z 6 x 3 i 



480 

z 3 x 5y2 

540 

x 3 zy 6 



+ 



+ 



+ 



540 

z 5 x 4 y 
480 

z 5 x 3y2 

540 

z 3 xy 6 



+ 



+ 



+ 



540 

z 5 xy 4 
480 

z 5 x 2y3 

540 

x 6 zy 3 



+ 



+ 



+ 



540 

x 4 zy 5 
480 

z 3 x 2yb 

540 

z 3 x 6 y 



+ 



xzy" 
540 
z 4 xy 5 
480 
z 2 x 3 y 5 
540 



z e xy 3 



30 

z 6 xy 



)N 



(III.1.40) 



+ 



+ 
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560 560 560 560 560 560 



z 7 xy 2 x 2 zy 7 z 2 xy 7 z 2 x 7 y z 7 x 2 y x 7 zy 2 
480 480 480 560 560 560 

A A 9 T A A O 



z 4 x 2 y 4 z 4 x 4 y 2 z 2 x 4 y 4 z 8 xy xzy 8 x 8 zy 
468 468 468 480 480 480 



..3„3..4 .,3„4,..3 _.9_.9..fi .,fi_9...9 ..9_.fi. .9)"^ 



2 



z 4 x 3 y 3 z 3 x 3 y 4 z 3 x 4 y 3 z 2 x 2 y & z & x 2 y 2 z 2 x®y 2 ' 

6300 6300 6300 6300 6300 6300 6300 6300 6300 

z 3 x 7 y 2 z 7 x 3 y 2 z 2 x 7 y 3 z 7 x 2 y 3 z 3 x 2 y 7 z 2 x 3 y 7 x 10 zy z 10 xy xzy 10 

6000 6000 6000 6000 6000 6000 6000 6000 6000 



n a a a o a a a n o a a ~t~ £ /i o ~t~ i k k £ o £ "I - 



z 4 x 2y6 z 2 x QyA z®x 2 y 4 z 4 x & y 2 z 2 x 4 y & z®x 4 y 2 z 2 x 5 y 5 z 5 x 2 y 5 z 5 x 5 y 2 
6300 6300 6300 6300 6300 6300 6660 6660 6660 



z 2 x 9 y x 9 zy 2 x 2 zy 9 z 9 xy 2 z 2 xy 9 z 9 x 2 y z®x 3 y 3 z 3 x 3 y® z 3 x®y 3 
6720 6720 6720 6720 6720 6720 5600 5600 5600 



~t~f7 /1~t~'7 A A *7 *7 A A H ~t~ OOO^OOO^ 



z 4 x 7 y x 7 zy 4 z 7 xy 4 z 4 xy 7 z 7 x 4 y x 4 zy 7 z 8 x 2 y 2 z 2 x 2 y 8 z 2 x 8 y 2 

6300 6300 6300 6300 6300 6300 | 5600 5600 5600 

z 2 x 9 y x 9 zy 2 x 2 zy 9 z 9 xy 2 z 2 xy 9 z 9 x 2 y z 8 x 2 y 2 z 2 x 2 y 8 z 2 x 8 y 2 

6720 6720 6720 6720 6720 6720 



z 8 x 3 y x 3 zy 8 x 8 zy 3 z 3 xy 8 z 8 xy 3 z 3 x 8 y 
6600 6600 6600 6600 6600 6600 



z 5 x 6 y z 6 xy 5 z 6 x 5 y z 5 xy Q x 6 zy 5 x 5 zy e 
6480 6480 6480 6480 6480 6480 6336 3 

z^x^y z^x^y^ z^x°y^ z^x^y z^x^y A z^x^y^ z^x^y* 



4 ) 

(III.1 



• Four-point functions: 

( |4), _ , 8 , 8 | 8 8 ; 8 8 

Pg y x i Vi z i w )~\'2 1^1 ? 22 2 2^ 2 2^22 

w^xy^z x^ywz* w L x £ yz xy^wz* w^xyz* x^y^wz 

6 6 6 6 

xy^wz xywz 6 w^xyz yx^wz 

72 72 72 ; 72 72 | 72 

w 3 x 2 yz 2 w 2 xy 3 z 2 w 2 xy 2 z 3 w 3 xy 2 z 2 w 2 x 3 y 2 z x 3 y 2 wz 2 

72 72 72 72 72 72 



w 2 x 3 yz 2 w 2 x 2 y 3 z x 2 y 3 wz 2 w 3 x 2 y 2 z x 2 y 2 wz 3 w 2 x 2 yz 3 
72 72 72 72 72 72 



x 2 ywz 4 


+ 


xy 4 wz 2 


+ 


w 2 xy 4 z 


+ 


w 4 xyz 2 


+ 


w 2 xyz 4 


+ 


w 2 x 4 yz 


72 


+ 


72 


+ 


72 


+ 


72 


+ 


72 


+ 


72 


1 4 9 
w^xy z z 


xy 2 wz 4 


x 4 y 2 wz 


x 2 y 4 wz 


w A x 2 yz 


x 4 ywz 2 


72 

+ ^ s- 


+ 


72 


+ 


72 


+ 


72 


+ 


72 


+ 


72 



w xy°z xy^wz w°xyz° x°y°wz w°x°yz x^ywz 

60 60 60 60 48 2 

yx°wz xywz w°xyz xy°wz w z x z y z z z 

576 576 576 ^ 576 576 576 



' w 3 x 3 y 2 z 2 w 2 x 3 y 3 z 2 w 3 x 2 y 3 z 2 w 2 x 2 y 3 z 3 w 2 x 3 y 2 z 3 w 3 x 2 y 2 z 3 
576 576 576 576 576 576 576 576 



w 3 x 2 yz 4 w A x 3 y 2 z x 3 y 2 wz 4 x 4 y 2 wz 3 w 3 x 4 y 2 z w 2 x 4 y 3 z x A y 3 wz 2 w 4 xy 2 z 3 
576 576 576 576 576 576 576 576 



~1~ a q 0~^~ 9 /)Q _ t~9/t Q Q A 9 ~t~ Q /I 9 9/1 Q ~t~ 9 Q /I 



w 4 x 3 yz 2 w 2 xy 4 z 3 x 2 y 4 wz 3 w 3 xy 4 z 2 w 3 x 4 yz 2 w 2 x 4 yz 3 w 2 xy 3 z 4 w 4 xy 3 z 2 



(III.1 
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576 576 576 576 576 576 576 576 

w 3 xy 2 z 4 w 4 x 2 y 3 z w 4 x 2 yz 3 w 2 x 3 yz 4 w 2 x 3 y 4 z x 3 y 4 wz 2 w 3 x 2 y 4 z x 2 y 3 wz 4 
480 480 480 480 480 480 480 480 

+ n r + — n K + —S 5- + -5 K + —a K + + —H a + 



w 2 x 6 yz xy 6 wz 2 w 6 xy 2 z ufixyz 2 x 6 ywz 2 x 6 y 2 wz w 2 xyz 6 xy 2 wz 6 
480 480 480 480 480 480 480 480 



w 6 x 2 yz x 2 ywz® w 2 xy 6 z x 2 y®wz w 5 xy 2 z 2 ufix 2 yz 2 x 2 y 2 wz 5 x 2 y 5 wz 2 

480 480 480 480 480 480 480 480 

w 2 x 2 yz 5 w 2 x 5 yz 2 ur 3 x 2 y 2 z w 2 x 2 y 5 z w 2 x 5 y 2 z x 5 y 2 wz 2 w 2 xy 5 z 2 w 2 xy 2 z 5 

540 540 540 540 540 540 540 540 



x 3 ywz 5 xy 5 wz 3 w 5 x 3 yz w 5 xyz 3 w 3 xy 5 z x 5 ywz 3 x 5 y 3 wz w 3 x 5 yz 

540 540 540 540 ^ 420 420 j 420 420 

w 5 xy 3 z x 3 y 5 wz xy 3 wz 5 w 3 xyz 5 xy 7 wz yx 7 wz xywz 7 w 7 xyz 



648 648 648 648 384 384 384 384 

w 3 x 3 y 3 z w 3 xy 3 z 3 x 3 y 3 wz 3 w 3 x 3 yz 3 w 2 x 2 y 2 z 4 w 4 x 2 y 2 z 2 w 2 x 4 y 2 z 2 w 2 x 2 y 4 z 2 
576 576 576 576 576 576 



x^y^wz w^x^yz x^ywz^ w^xy^z w^xyz^ xy^wz^ 



1.6 Recurrent formulas for Gaussian correlators 

• General recurrent formula for connected correlators with tr^>: 13 



(III.1.43) 



G 

ai!(ti-oi)! a m \(i m -a m )\ W r Hg 



0<ai ,...,a m <fc 
ai+...+a m =fc 



Particular examples: 



(W>) fc )) G = ,5 fc ,o + iV5 fc , 2 (HI.1.45) 



(tr^) K tTtf)) a = N z 5 kfi + 2N5 K2 (III.1.46) 
(tr# tr^ 3 )) G = 3iV 2 5 fcil + 6N8 k , 3 (III.1.47) 



(tr0) K trc/) 4 )^ = (2iV d + iV)<S fc)0 + 12JV^ fc)2 + 24/V<S M (III.1.48) 

(((tr</>) fe tr0 fc )) c = ATA;! (III.1.49) 

(((W>) fc tr^ m )) G = 0, for k > m (III.1.50) 
In general, for the k + 1-point correlators 

; (tr# tr^ 2 ™))^ = c 2m <5 M + 2m(2m - l)c 2m _ 2 <5 fe , 2 + . . . + (2m)\N5 k , 2m (III.1.51) 



(tr0)* tr^ m+i ^ c = (2m + l)c 2n A,i + 2m(4m" - l)c 2m _ 2 <5 fe)3 + . . . + (2m + l)!iV<5 fe , 2m+ i 



(III.1.52) 



In applications of this formula one should be careful: when SOH16 %a. &a = 0, while some other ip — ap > 0, the 
corresponding connected correlator at the r.h.s. . . . ((tr/ . . . tr0 l < 3 ~ a ' 3 ■ ■ -)) G = 0- 
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where 



C2m = ((trcj) 



,,2/v 



The disconnected correlators 



G 



(ti^f k ) = (2k - l)UN k 

I G 



• Gaussian correlators of tr0 2 only 



(W> 2 )^ = 4(fe - l)(fc - 2) ^(tr0 2 ) fc ^ = 2 fc ~ 1 (A; — 1)!AT 2 



(tr</> 2 ) 



• 2(/c - 1) \N 2 + 2(A; - 2)] ^(W 2 ) 



k-2 



G 



E^r E 



A:! 



— j! — 1 X \...IA 



2 VJ 



G 



h+...+L=k 



E^ 2i E 



Particular examples: 



0<! 1 ,...,!j<*: 
h + ...+«j=fc 



(tr0 2 ) c = N 2 



u • ■ ■ ■ • u 



(W 2 ) : 



iV 4 + 2iV 2 



G 



(W 2 )' 



iV 6 + 6N 4 + 8iV 2 



G 



tr^J ^ = N a + 12iV° + 44iV 4 + A8N Z 
(trc/> 2 ) 5 ^ = iV 10 + 20iV 8 + 140iV 6 + 40(W 4 + 384iV 2 



(III.1.53) 



(III.1.54) 



(III.1.55) 



(ni.1.56) 
(III.1.57) 

(III.1.58) 

(III.1.59) 

(III.1.60) 



• Recurrent formulas for correlators with ti(f> k with k > 2 

These are more sophisticated, from the point of view of iV dependence (genus expansion). See 
pil.l.lOfiJ) and (|III.1.120J> for particular examples with tr0 2 . 

1.7 Lowest Gaussian correlators, connected and disconnected 

To restore the g and M dependencies below, each correlator 

(ti4> kl . . . tT(p km ^ 

should be multiplied by 



r 1 + ... + K m 



M . 



• Selection rules: 

a) for disconnected correlators for N = 1 



tr^ 1 ...tr^) G = (h + ...+i k -l)\\ 



(III.1.61) 



b) i\ + . . . + ik = even 
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Single-point correlators (the generalized Catalan numbers): 



^ 2 ) G = (K)) G = Ar2 



W> = ^ //fi 



5iV 4 + ION 2 



tT(ff) G = {{tup}} = UN + 7(W d + 217V 
tr0 lo ) c = ((tr0 lo )) G = 42N 6 + 420iV 4 + 483iV 2 



Two-point correlators: 



(tr^) 2 )) G = ((tr^) 2 ) G = N 



I rrr ) )) = 2N 2 



(trcff 



N 4 + 2N 2 



G 



(tr</> 3 ) \\ = /( tr0 3 ) \ = 12iV 3 + 3iV 



(tr^ 4 ) : 



36iV 4 + 6(W 2 



G 



(tr^ 4 ) * ) = 4iV 6 + 40iV 4 + 61iV 2 



G 



(tr^ 5 ) : 



G 



(tr0 5 ) ^ = 1807V 5 + 6007V 3 + 1657V 



tr0 tr^ ) ) = (tnf) to(f> ) = 37V^ 



/.5\ 



tr0 trcj) } = 107V' 4 + 57V 



ti(f) te<P 7 )) G = (W> te<P 7 ) G = 357V 4 + 707V 2 



G 



tr^> W> 9 ))^ = (tr</> tr^ 9 ) = 1267V 5 + 6307V 3 + 1897V 



tr0 2 tr0 4 )) = 8 A*'"' + I V 



tr0 2 tr0 4 ^> G = 2iV 5 + 97V 3 + 4iV 



tr^ 2 tr</> 6 



G 



30iV 4 + 60iV 2 



tr</> 2 tr0 6 ) = 57V b + 407V 4 + 607V 2 

G 



tr^ 2 tr^ 8 ^ = 112iV 5 + 560iV 3 + 1787V 



tr0 2 tr# 8 ^ = 14iV 7 + 1827V 5 + 5817V 3 + 1787V 
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trcj) 3 tr^ 5 ))^ = (tr0 3 tr0 5 ) G = 45iV 4 + QON 



tr</> 3 ^<P 7 )) G = (trcp 3 tr<p 7 ) G = 168iV 5 + 63(W 3 + 147iV 



tr0 4 tr</» 6 ^ = 144iV 5 + 600iV 3 + 156iV 

tr<p 4 tr<p e ) = WN 7 + 169iV 5 + 61(W 3 + 1567V 

/ G 



Three-point correlators: 



(w> 2 ) ; 



8N 2 



G 



(tr^ 2 )' 



N 6 + 6iV 4 + 87V 2 



G 



(W>) 2 tr^ 2 )) G =2.V 



(tr</>) 2 tr0 2 ) =iV 3 + 2iV 



(tr</>) 2 W* 4 ^ = 12iV 2 
((W>) 2 tr</. 4 ) G = 2iV 4 + 13iV 2 



(trcj)) 2 trtf )) - 00.V' + 'InV 



(tr0) 2 tr/) = 5iV 5 + 7(W 3 + 30N 
I G 



(tr^) 2 te<f)) G = 280iV 4 + 560iV 2 



{ti(j)) 2 tr/\ = 14iV 6 + 350iV 4 + 581iV 2 



(tr</> 2 ) W> 4 ^ = 48iV 3 + 24iV 
(W 2 ^ tr^ 4 \ = 2N 7 + 21iV 5 + 58iV 3 + 24iV 



6 \ \ = 24(W 4 + 4807V 2 
G 



(tr^ 2 ) tr/^ = 5N 8 + 8(W 6 + 38(W 4 + 480iV 2 



(tr</> 3 ) W> 2 ^ = 72iV 3 + 18iV 
(tr^ 3 ) 2 tr0 2 ^ = 12iV 5 + 75iV 3 + 187V 
(tr0 3 ) 2 tr^ 4 ^ = 432iV 4 + 468iV 2 
tr</> 3 ) 2 tr0 4 \ = 24iV 6 + 45(W 4 + 471iV 2 



(III.1.84) 
(III.1.85) 

(III.1.86) 
(III.1.87) 

(III.1.88) 
(III.1.89) 

(III.1.90) 
(III.1.91) 
(III.1.92) 
(III.1.93) 
(III.1.94) 
(III.1.95) 
(III.1.96) 
(III.1.97) 

(III.1.98) 
(III.1.99) 
(III.1.100) 
(III.1.101) 

(III.1.102) 
(III.1.103) 
(III.1.104) 
(III.1.105) 
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(tr</> 4 ) ti(f> 2 

2(8iv((tr0 2 tr0 4 )) + 4 (((W>) 2 W> 4 )) ) + 2 • 16 ((W> 6 )) = (III.1.106) 



G W v I I G' W 1 1 G 

288iV 4 + 480iV 2 



(tr0 4 ) tr0 2 ^ = 36iV 6 + 36(W 4 + 549iV 2 (III.1.107) 

^tr0 tr^ 2 tr^ 3 ^ G = 12iV 2 (III.1.108) 

tr</» tr</> 2 tr0 3 \ = 3iV 4 + 12iV 2 (III.1.109) 

/ G 

tY<t> tY<t> 2 tr</> 5 )) G = 60iV 3 + SON (III.1.110) 

tr0 tr0 2 tr^> 5 \ = 10iV 5 + 65iV 3 + 30 N (III.l.lll) 



G 

tr0 tr</> 2 tr^ 7 ^^ = 28(W 4 + 56(W 2 (III.1.112) 

tr</» tr^ 2 tr</> 7 ) = 35iV 6 + 350iV 4 + 560iV 2 (III.1.113) 

/ G 



tr0 tr0 3 tr(/> 4 )) G = 72iV 3 + 24iV (III.1.114) 



tr<p tr^ 3 tr</> 4 ) = 6iV 5 + 75iV 3 + 24iV (III. 1.115) 

tr0 tr</> 3 tr0 6 )) G = 36(W 4 + 540N 2 (III.1.116) 

tr<£ tr0 3 tr<£ 6 ) = 15iV 6 + 390iV 4 + 540iV 2 (III.1.117) 

/ G 

(Jjx<j> tr</> 4 tr0 5 )) G = 36(W 4 + 540N 2 (III.1.118) 

(iicp tr0 4 tr^ 5 ) G = 2(W 6 + 380iV 4 + 545iV 2 (III.1.119) 

tr^ 2 tr^ 3 ^)) G = 6N ((trcf) tr</> 5 )) G + 

+ (\(W^(tre/> 3 ) 2 ^ + io((tr0 tr^ 2 trc/. 3 ))^ + 2 ■ 3 ■ 5 ((tr^ 6 )) G = (III.1.120) 

= 36(W 4 + 480iV 2 

(trcj) 2 tr0 3 te(f) G = 45iV 6 + 420iV 4 + 480iV 2 (III.1.121) 
Four-point correlators: 

\ti4)f)) = (III.1.122) 



G 

{tr<j)) A ) = 3N 2 (III.1.123) 



G 

4 



(tr0 2 ) J J =48N 2 (III. 1.124) 
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(W 2 y 



TV 8 + 12/V 6 + 447V 4 + 48iV 2 



G 



(((W>) 3 tr^ 3 )) G = 6iV 



(tr0) 3 tr^ 3 ^ = 9N 3 + 6N 
(((tr^) 3 tr^ 5 )) G = 60iV 2 
(tr^) 3 tr^ 5 ) G = 30iV 4 + 75iV 2 
(tr^) 3 te<P 7 )) G = 420iV 3 + 210iV 
((W>) 3 tr0 7 ) c = 105iV 5 + 630iV 3 + 210iV 

^(tr^. 2 ) 3 tr^ 4 ^ = 384/V 3 + 192iV 
(tr^ 2 ) 3 tv(f> 4 ^ = 2N 9 + 37iV 7 + 226iV 5 + 488iV 3 + 192iV 
(tr</> 3 ) 3 tr^ = 648iV 3 + 162iV 
(tr0 3 ) 3 tr<^ = 108iV 5 + 675iV 3 + 162iV 



8N 



G 



(ti4>) 2 (ti4> 



2 ^ \ = AT 5 + 6iV 3 + 8N 
G 

2\ 



(tr0) 2 (tr</> 3 ) \ \ = 72iV 



/(W>) 2 (tr<j? 



3(W 4 + 75N' 



G 



(tr0) 2 (tr</> 4 ) J J = 576iV 3 + 192iV 



(tr0) 2 (tr</> 4 ) 2 ^ = AN 1 + 88N 5 + 661iV 3 + 192iV 



(tr0 2 ) (tr<^ 3 ) 
(tr(p 2 y (trcf) 



576iV 3 + 144iV 



G 



3 ^ 2 \ - 12iV 7 + 171iV 5 + 618iV 3 + 144iV 



{tY<t>Y tr<^ tr^ 4 ))^ = 72A^ 



(tr^) 2 tr^> 2 tr0 4 ) = 2N 6 + 257V 4 + 787V 2 



(tr<f>) 2 trcf> 2 tr/)) - 180 V ; + 2 10.Y 



ni.i.125) 
ni.i.126) 

QI.1.127) 
QI.1.128) 
QI.1.129) 
III.1.130) 
III.1.131) 

[II. 1.132) 
QI.1.133) 
HI.1.134) 
QI.1.135) 

QI.1.136) 
QI.1.137) 
QI.1.138) 
III.1.139) 
QI.1.140) 
QI.1.141) 

QI.1.142) 
QI.1.143) 

[II. 1.144) 
QI.1.145) 
HI. 1.146) 
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(tr</>) 2 tret? ticf) = 5N 7 + 11(W 5 + 59(W 3 + 240JV 

I G 



tc<t> (tr^ 2 ) ticfjj = 72N 2 
ti<t> (tr</> 2 ) 2 ti(f> 3 \ = 3iV 6 + 30iV 4 + 72N 2 



G 

tT(f> (tr^ 2 ) 2 W> 5 ^ = 480iV 3 + 240iV 
tr(f> (tr0 2 ) 2 tr<f> 5 ^ = WN 7 + 145iV 5 + 55(W 3 + 240N 



(tr</>) 2 tr(j> 3 tnfP)) = ■> IDA"' + 210A' 



(tr</>) 2 trcj) 3 tr^ 5 ) = 105iV 5 + 630iV 3 + 2107V 

/ G 



ti(f> trcj) 2 ti<j) 3 tr</> 4 )) G = 576iV 3 + 192iV 



ticf) ti(j) 2 trip 3 tr^> 4 ^ G = 126iV 5 + 627iV 3 + 192iV 



• Five-point correlators: 



384iV 2 



G 



(W 2 )'' 



N 10 + 20iV 8 + 140iV 6 + 40(W 4 + 3847V 2 



G 



(tr^) 4 tr0 2 )) G = 



(tr0) 4 tr0 2 
(tr(/>) 4 ti(p 



G 

4 4-„^4 



3iV 4 + 12iV 2 



G 



24iV 



(tr^) 4 W* 4 ^ = 6iV 5 + 75iV 3 + 24iV 



(tr</>) tr<^ tr^)) = 36N 



(ti(j)) 3 tr<f> 2 tr<f> 3 ) = 9N 5 + 60iV 3 + 36iV 

/ G 



(tr^) 3 tr0 2 tr^ 5 ^ G = 48(W 2 
((W>) 3 tr</> 2 tr0 5 ) G = 30iV 6 + 315iV 4 + 60(W 2 



(tr0) 3 tr0 3 tr^ 4 ^ G = 504iV 2 



(tr0) tr^ tr0 4 ) = 18iV° + 3457V 4 + 582iV z 

/ G 
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I it;)- ) ) = 48N 



(tr^ 2 ) (tr 

(tr0 2 ) d (tr(f>) 2 \ = N 7 + 12N 5 + UN 3 + A8N 



(tr<f>) 2 (tr^ 2 ) 2 W> 4 \\ = 384iV 2 



(tr0) 2 (W 2 ) 2 tr^A = 2N 8 + 41iV 6 + 278iV 4 + 624iV 2 



(III.1.168) 
(III. 1.169) 

(III.1.170) 
(III.1.171) 



Six-point correlators: 



(tr0) D ) G = 15iV d 



(III.1.172) 
(III.1.173) 



(tr0) 5 tr0 3 )) c = 



(tr^) 5 tr^ 3 ) G = 45iV 4 + 



60N 



^5 +T ,^3\ 

(tr0) 5 tr</> 5 )) G = 120 A' 



(tr0) 5 W* 5 ^ = 150iV 5 + 675iV 3 + 120iV 



(tr0) 4 (tr0 2 ) )) =0 



(tr</>) 4 (t^ 2 )^ 



3iV 6 + 30iV 4 + 72iV 2 



(tr^) 4 ftr^ 3 "' 



= 216iV 



G 



(tr^) 4 (tr0 3 ) 2 ^ = 144iV 5 + 585iV 3 + 216iV 



(tr</> 2 ) 4 (tr(py 



384iV 



G 



(tr</> 2 ) 4 (tr^) 2 \ = iV 9 + 20iV 7 + 140iV 5 + 40(W 3 + 384iV 



(tr^) 3 (trc^y tr^ 



288iV 



G 



(tr</>) 3 (tr^ 2 ) tr^y = 9N 7 + 132iV 5 + 516iV 3 + 288iV 



(III.1.174) 
(III.1.175) 
(III.1.176) 
(III.1.177) 

(III.1.178) 
(III.1.179) 

(III.1.180) 
(III.1.181) 

(III.1.182) 
(III.1.183) 

(III.1.184) 
(III.1.185) 



• Seven-point correlators: 



(tr^) 6 tr^ 2 )) G = 



(tr^) 6 tr^ 2 ) G = 15iV° + 



90iV 3 



(III.1.186) 
(III.1.187) 
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(tr^) 6 W> 4 )) g = 



(trc/)) 6 tr^ 4 ) = 30iV 6 + 555iV 4 + 360iV 2 

/ G 



(W>) 5 tr^ 2 tr0 3 )) G = 
(tr0) 5 tr</> 2 tr0 3 ) = 45iV 6 + 42(W 4 + 48(W 2 

/ G 



(W>) 4 (W 2 ) ^ =0 



(tr0) 4 (tr^> 2 ) 3 ^ = 3N 8 + 54iV 6 + 312iV 4 + 576iV 2 



Eight-point correlators: 



(t^T)) G = o 



(tr</>) 8 ) G = 105iV 



(tr^) 7 tr^)) G = 



(tr^) 7 tr^ 3 ^ = 315iV 5 + 63(W 3 



(tr0) 6 (tr^ 2 ) 1 



G 



(tr0) 6 (tr</> 2 ) ^ = 15iV 7 + 210iV 5 + 720iV 3 



Nine-point correlators: 



(W>) 8 tr^ 2 )) G = 



(tr^) 8 tr0 2 ^ G = 105iV 6 + 840iV 4 



Ten-point correlators: 



(tr0) lo ) G = 945iV 



47 



1.8 Lowest prepotentials J 7 ^: first terms of i-expansions 

The prepotential is the generating function of the connected correlators 

oo -.00 

^ G = A^ o + £t fc ((tr0 fc )) +- £ t fcl t fc2 ((tr^tr^)\ + 



G 2! ^ \\ r / / G 

fc=0 fci,fc 2 =0 

^ (III. 1.204) 

+ £j E **i *fa *fc 3 ( (tr</> fcl tr<^ 2 tr0 fc3 ) ) q + . . . 

ki,k2,k3=0 

If the correlators are taken from Table IT771 one should restore g and M dependencies and substitute 

N = — (III.1.205) 

9 



so that 



F G {t\S) = Y J 9 2p - 2 ^\t\S) 

OO OO fc]_ +• •■ + fcm 



, ml , V „ <? m \M J W //G\ T IG 

m=l ki,...,km=0 a 



Nt N 2 t 2 g(2N 3 + N)t 4 g 2 (5N 4 + WN 2 )t 6 g 3 (UN 5 + 70N 3 + 21N)t 8 
~g W + M 2 + M 3 + M 4 + 

g 4 (42N 6 + 42(W 4 + 4837V 2 )t 10 
+ W> + " " " + 



Nt 2 N 2 t 2 2 g(12N 3 + 3N)tl # 2 (36iV 4 + 60N 2 )t 2 4 g 3 {180N 5 + 600iV 3 + 165iV)t§ 



+ -rrrf + — t"t~7o ^ + ^ — — , + — :_£ + ...+ 



2 5 M M 2 2M 3 2M 4 2M 5 



3iV 2 tit 3 ff (10iV 3 + 5N)tit s 5 2 (35iV 4 + 70A^ 2 )tit 7 5 3 (126A^ 5 + 6307V 3 + 189N)ht 9 
+ M 2 + M 3 + M 4 + W> + " ' + 

, g{8N 3 + AN)t 2 U + ff 2 (30iV 4 + 60iV 2 )t 2 t 6 + ff 3 (112iV 5 + 560iV 3 + 178N)t 2 t 8 + 



M 3 M 4 M 5 

5 2 (45A^ 4 + 60iV 2 )t3t5 5 3 (168iV 5 + 630A^ 3 + U7N)t 3 t 7 
+ M 4 + M5 + -- - + 

5 3 (144iV 5 + 6007V 3 + 156iV)t 4 t 6 8iV 2 t| 
+ M* + " ' + ~6M 3 ~ + " ' + 

2Nt 2 t 2 l2N 2 t\t A g(60N 3 + 30N)tjt 6 g 2 (280N 4 + 560N 2 )tjt 8 
' 2gM 2 H 2M 3 + 2M 1 + + • • • + 

g(A8N 3 + 24iV)t|t 4 5 2 (240iV 4 + 480iV 2 )i 2 i 6 
+ 2M 4 + 2M~ 5 + ■ • • + 

<?(72iV 3 + 18A%i 2 5 2 (432iV 4 + 468A^ 2 )t§t 4 
+ 2M 4 + 2M5 + ■ • • + 

5 2 (288iV 4 + 480iV 2 )t 2 i4 
+ 2M5 + '" + 

lgjgfifgfg g(60iV 3 + 30AQtit 2 t 5 g 2 (280iV 4 + 560iV 2 )M 2 t 7 
+ M 3 + M 4 + M5 + • • • + 

(III.1.206) 
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g(72N 3 + 24A0Ms*4 <? 2 (360iV 4 + 540N 2 )ht 3 t e 

+ w + w + ■ • • + 

g 2 (36(W 4 + 54QjV 2 )tit 4 t5 
+ IP + • • • + 

g 2 (360iV 4 + 480iV 2 )t 2 t3t5 , fSiV 2 ^ 

+ M 5 + '"" + 24M 4 + "' + 

6iVtft 3 60AT 2 t3t 5 5 (42(W 3 + 2WN)t\t 7 <?(384iV 3 + 192iV)^ 4 

+ 63M 3 + 6M 4 + 6M5 + • • • + ^5 + • • • + 

SNtjtl 72N 2 tjtl g(576N 3 + !92N)t 2 tj 
+ AgM 3 + 4M 4 + 4M5 + • • • + 

g(5767V 3 + lUN)t 2 t 2 
+ 4M5 + ' • • + 

72N 2 t\t 2 t± g(540N 3 + 210N)tlt 3 t 5 5 (480iV 3 + 240N)t 2 t 2 t 6 
+ 2M 3 + 2M5 h 2M5 + • • • + 

72iY 2 M 2 ,t 3 g(48(W 3 + 240iV)M 2 t 5 
+ 2M 4 + 2M5 + ••• + 

g(576iV 3 + 192jV)tit 2 *3*4 
4 M5 ■ ■ ■ + 



384iv2 ,5 (HI.1.207) 

H + + 

120M 5 

24iVi 4 t 4 



24M 5 



■ + ...+ 



36Ntft 2 t 3 480N 2 tlt 2 t 5 504iV 2 t 3 t 3 t 4 
+ 65M 4 + 6M5 + 6M5 +••• + 

48Nt 2 4 
H — + + 

384iV 2 £ 2 i 2 ,i 4 

; -+...+ 



4M 5 
-I LJi _i_ _i_ 



216iVt 4 t| 
+ 48#M 5 

384iV£|i| 
+ 48#M 5 

288Nt 3 t 2 .t 3 

H + . . . + 

I65M 5 



+ ...+ 
+ ...+ 



49 



1.9 Functions Za(t\N): first terms of t-expansion 

e- Nt »Z G (t\N) = 1+ 

N 2 t 2 g(2N 3 + N)t 4 g 2 (5N 4 + WN 2 )t 6 g 3 (UN 5 + 70iV 3 + 21N)t 8 
+ ~M~ + M 2 + M 3 + M 4 + 

g 4 (A2N 6 + 420iV 4 + 483iV 2 )ti 

+ w> + • • • + 

Nt\ (N 4 + 2N 2 )t 2 2 g{\2N 3 + 37V)t 3 g 2 (4N (i + AON 4 + 6lN 2 )tl t 



2gM 2M 2 2M 3 2M 4 

g 3 (180N 5 + 600iV 3 + 165iV)t§ 



2M 5 



+ ...+ 



3jV 2 tit 3 g(10iV 3 + 5N)ht 5 g 2 (35N 4 + 70iV 2 )tit 7 
+ M 2 + M 3 + M 4 + 

g 3 (126iV 5 + 630iV 3 + 189AQM 9 

+ W> + "" + 

g(2N 5 + 9N 3 + 4A%i 4 g 2 (5N 6 + 40iV 4 + 60N 2 )t 2 t 6 
+ W + W 4 + 

g 3 (UN 7 + 182iV 5 + 581iV 3 + 178N)t 2 t 8 
+ W> + "" + 

g 2 (45iV 4 + 60iV 2 )t3t 5 ff 3 (168iV 5 + 630iV 3 + 147AQt 3 t 7 
+ M 4 + W> + " ' + 

g 3 (WN 7 + 169iV 5 + 610iV 3 + 156N)t A t 6 
+ W> + '"" + 

(iV 6 + 6iV 4 + 8N 2 )4 
+ 6M 3 + "' + 

| (AT 3 + 2N)tjt 2 (2N 4 + l3N 2 )t 2 t 4 g{5N 5 + 70iV 3 + 307V)*? t 6 



2gM 2 2M 3 2M 4 

g 2 (UN 6 + 350iV 4 + 5817V 2 )i 2 i 8 
+ 2M5 + • • • + 

g(2N 7 + 21iV 5 + 58iV 3 + 2AN)t 2 2 U g 2 (5N 8 + 80iV 6 + 380iV 4 + 480iV 2 )t|i 6 ^ 



2.1/ 1 2M 5 
g(12N 5 + 75N 3 + l8N)t 2 tl g 2 (2AN 6 + 450iV 4 + 4717V 2 )*§t 4 
+ 2M 4 + 2M5 + • • • + 

g 2 (36iV 6 + 360iV 4 + 5497V 2 )t 2 tj 
+ 2M5 + • • • + 

(3iV 4 + 12iy 2 )tit 2 t 3 g (10iV 5 + 657V 3 + 3QjV)tit2*5 
+ M 3 + M 4 + 

<? 2 (35iV 6 + 350iV 4 + 560iV 2 )tit 2 t7 

+ w> + ■ • • + 

, g{6N 5 + 75iV 3 + 24iV)iit 3 t 4 g 2 (15iV 6 + 390iV 4 + 540iV 2 )M 3 t 6 ^ 



M 4 M 5 
g 2 (20N 6 + 380iV 4 + 545jV 2 )tit 4 t 5 
+ IP + ■ • • + 

<? 2 (45iV 6 + 420iV 4 + A80N 2 )t 2 t 3 t 5 

+ w> + ■ • • + 
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3N 2 ti N s + UN 6 + UN 4 + 48N 2 )d 

-I ± 1 i-^. + + 

24g 2 M 2 24M 4 -r . . . -i- 

(9iV 3 + 6iV)tf i 3 (30iV 4 + 75N 2 )tft 5 g(W5N 5 + 630iV 3 + 21(W)tft 7 
+ %AP + 6M 4 + 6M5 +••• + 

c/(2iV 9 + 37N 7 + 226iV 5 + 488iV 3 + 192N)4U 
+ 6M5 + • • • + 

(iV 5 + 6iV 3 + 8N)t\t 2 2 (30iV 4 + 75iV 2 )t 2 t§ g{AN 7 + 88iV 5 + 661iV 3 + 192iV)t?^ 



+ - 7TT2 '~^ + - : — = : ^^ + ...+ 



4c/M 3 4M 4 4M 5 

#(12iV 7 + 171iV 5 + 618iV 3 + U4N)44 
+ 4M5 + ■ ■ ■ + 

(2iV 6 + 25iV 4 + 78N 2 )t\t 2 U 



+- 



2M 4 



g(105N 5 + 630iV 3 + 210N)tjt 3 t 5 g{5N 7 + HOiV 5 + 59(W 3 + 240N)tj t 2 t 6 
+ 2M5 + 2AP + • • • + 

, (3iV 6 + 3(W 4 + 72iV 2 )M^3 g(l(W 7 + 145iV 5 + 55(W 3 + 240AQt 1 t 2 .t 5 



2M 4 2M 5 
g(l26N 5 + 627iV 3 + 192^)^2*3*4 

+ w> + ■ • ■ + 

(iV 10 + 20iV 8 + UON 6 + 40(W 4 + 384iV 2 )^ 

+- — + + 

120M 5 

(3iV 4 + YlN 2 )t\t 2 (6iV 5 + 75iV 3 + 2AN)tjt 4 
+ 24g~M* h 24M 5 +••• + 

(9iV 5 + QON 3 + 36N)tft 2 t 3 (30iV 6 + 315iV 4 + 600N 2 )tft 2 t 5 
+ 6gM 4 + 6M5 + 

(18iV 6 + 345iV 4 + 582N 2 )tft 3 t 4 
+ 6M5 + ' • • + 

(TV 7 + 12iV 5 + 44iV 3 + A8N)tl4 
+ \2g~M 1 +, " + 

(2iV 8 + 41iV 6 + 278iV 4 + 6247V 2 



AM 5 
15N 3 4 
+ 720<? 3 M 3 + "' + 

(457V 4 + 60AQtf t 3 (150iV 5 + 6757V 3 + 120jV)tfo 5 
+ 120c/ 2 M 4 + I2O5M 5 +••• + 

(3/V 6 + 3(W 4 + 72N 2 )t\t 2 2 (144iV 5 + 585iV 3 + 216iV)tf*§ 



+ - T~ — TTk + 



48# 2 M 4 48#M 5 

(iV 9 + 20iV 7 + UON 5 + 40(W 3 + 3847V)t 2 t| 
+ 48g~AP + • • • + 

(9iV 7 + 132iV 5 + 5167V 3 + 288jV)tfo|t 3 
+ 16gM 5 + • • • + 

(15iV 5 + 90N 3 )4t 2 + (30iV 6 + 555iV 4 + 360iV 2 )^t 4 + + 



720# 3 M 4 720# 2 M 5 

(45iV 6 + 420iV 4 + A80N 2 )tlt 2 t 3 
+ 60 5 2 M 5 + • • • + 



(3iV 8 + 54iV 6 + 312iV 4 + 576N 2 )tf4 
240# 2 M 5 
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+ ...+ 



(III.1.209) 



105iV 4 if | (315AT 5 + 630N 3 )tlt 3 

+ 5760c? 4 M 4 + ' ' ' H 720c/ 3 M 5 + • • • + 

(15iV 7 + 210iV 5 + 720iV 3 )ift| 
+ 180 5 3 M 5 + " ' + 

(105iV 6 + 840N*)tft 2 (III.1.210) 
5760# 4 M 5 + • • • + 

945iV 5 t}° 

H — + + 

10!<? 5 Af 5 

+ ... 

The MAPLE generating program for the expansion of the integral 



j"#e-±Trc/> 2 



(III.1.211) 



in g and written by V.Pestun, is available at http://wwwth.itep.ru/progs/giv.mws or .../giv.tar. It 
generates the function "tttr", which is the integral (|III.1.211|) and requires, at least, MAPLE 7. 

2 Non-Gaussian partition functions 

From the Virasoro constraints in the case of non-Gaussian potential one can get non-Gaussian densities 
and multidensities, for instance, the simplest genus zero one-point density is as follows 

P T\z) = W ' {z) ~ v{z) On.2.1) 

Here 

y 2 w (z)=W'(z) 2 -lf$ 1 \z) (III.2.2) 
From equation (jl.3.40f) follows expressions for all non-Gaussian (multi) densities, for instance, 

(0|1), x (0|1), 



and 



= ^(p$V|*) (ni.2.4) 

It can be shown explicitly that the function p^) 2 \z, x) is symmetric in z and x, if /^J' 2 ^(z|a;) is defined 

as in ai.a.am 

/g |2) (z|x) = ii w (z) /0 SJ 1) (x) (III.2.5) 
The two point function can be presented in form explicitly symmetric in x and z, namely, 

pT\ z ^) = 9 7 w w ^(W(x)W(^) -yOr)y(z) + (x - Z f%( v )F (0) + (0) ) 

2y(x)y{z){x — z) z \ / 

(III.2.6) 

where 

n— 2 m— 2 ^ g 

K w (x,z)= £ nmTX^^sV- — (III.2.7) 

n,m=0 i=0 j=0 01 n -i-2 01 m -j-2 



and 



r) n— 2 a 

A w (x,z) = 2 $>T n (2— - + - i){x - z)(x^ 1 - z*" 1 ) + (zx 4 " 1 + xz 1 "' - ' 



(III.2.8) 
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Thus, the two-point function has a smooth limit when x — > z 

M2), 1 /I 



Pw {x,x) 



2y(x) 2 V8 



(Uw^W -2W" 2 ) 
\ 8 



YT^((( W ") 2 + WW® - 2R'{ v (z)F^)(W 2 - 4R' w (z)F^) - 2(W'W" - 2R' w {z)F^) 2 ) 

+K w {x,x)F^ +B(x)F^)) 

n-2 



1 



B(x) = o E nT n E 4 (" - *')(* - 2) + (* - 1)(* - 4 



j~2 







1=1 



ST, 



n— i— 2 



(III.2.9) 
(III.2.10) 



«trj»w=x;^ 



i=l 

tr<^ fc )) T j = J2 for < k < n, 



i=l 



(m.2.11) 



tr^+ 1 )y ( ° ) =f^5 i <+ 1 + ' 



i=i 



ra+l)T„+i ^ 



3 CIV-DV prepotential: Q through Gaussian prepotentials 

From 1)11. 2.5(1 we obtain here (as a generalization of the results of |23| ) the first terms of expansion 
for the prepotential through Gaussian prepotentials. They coincide with the results from 

Seiberg-Witten theory (see [571] . [BT] . ) : 



lt=o 



t=o 



^(S 2 - 45i5j + S]) log ay - ^ Sf+ 



i oo 2fc-2 

F 2 -(5a)F 1 + -^^ r (c95) fc F 1 2 
fc>i 



DV 



i=0 



J F 1 3-(^)F 2 + -(5a) 2 F 1 + 



oo „2fc-2 



fc>l 



ddf [F\F 2 - F 1 {Sd)F 1 ] + - V 9 -—- ( (cdd) k F? - ZF^cddfF 2 

6 — ' fc ^ 

fc>2 



-r 4 



t=0 



F A - (Sd)F 3 + \{Sd) 2 F 2 - \(SdfF x + 

2 D 



Ai)__rp(i) 



2fc-2 



fc>i 



V 2 + F X F Z - F 2 (59)F! - F l {Sd)F 2 + ^F 1 {Sd) 2 F 1 + X - \{Sd)F x f 



+ 



oo 2fc-4 



-1 LXJ 

1 5 



+ 2^ fcl 



{cdd) k F 2 F 2 - 2F 1 (cdd) k F 1 F 2 - F 2 (cdd) k F{ 



k rp2 



+ 



°° „2fc-4 r 



{cdd) k F 2 {Sd)F 1 - 2F 1 {cddf F l {Sd)F 1 - (Sd)F l (cdd) k F{ 



k jp2 



+ 



(III.3.1) 
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1 9 



-3 {{cdd) 



(cdd) k F? + 12F?(cdd) k F? - AF^cddfFf 

i ® (caa)) k Ft ® Ft 



+ 



(p|m) 



1 Explicit formulas for JF DV 
= 1 : 



t=o 
(o|i) 



DV 



t=0 



«4 1 )) (1) m , ((3 2 )) (1) 

a 2 + 2! . 32 




,(o 



3(7. 



(0 



^DV 



4"+ 



«3',7')) m ,„ <„ «4',6'))' 2 ' <„ (i) «5'» ( ' 2) 

— O J it ~ | •' ' ' 



4 3-7 

■2 «1\\(2) . 



4-6 



<w + 



2! -5 2 



(7 



(0 



2! • 3 2 • 6 

«3 2 ,4 2 ))< 2 > 
(2!) 2 • 3 2 • 4 2 



(0 



(0 «3 1 . 4l .5 1 )) (2) .(0.(0.(0 «4 3 » 



(2) 



(7 



(0 



2 r a g> 



3-4-5 



a 1 a 2 ct 3 



3! -4 3 



a 2 



+ 



2 ((S 4 ,^)) 



(2) 



4! • 3 4 • 4 



a 



(0 



(2) 



((3 6 )) 
6! -3 6 



(7 



(0 



m = 2 



(0|2) 



(1|2) 



t=0 



DV 



E 



= 4 E( 5 * - + S D lo § «v - ^ E 

f , ((4 2 )) (1) r (Q 



M 2 ^ 6 

((3 2 ,4)) (1) 
2! • 3 2 • 4 



(7 



(0 



2! -4 2 

(0 ■ ((3 4 )) (1) 



3-5 



1 u 3 



oV + 



4! • 3 4 



(7 



(0 



V S,S 3 - / ((4 1 )) (1) 2((3 2 )) (1) (<) «2,4))W w ((2!,3 2 )) (1) 
^2" I ?m 1 °l 771 "2 



24 



9a?- 



4a 2 - 



, 2«1,5))W (Q 2«1,3,4))W (i) (0 2«1,3 3 )) (1) 
+ 5a t , ^ + 3.4a,, ° x + 3! • 33a,,- 



2! • 3 2 a 2 - 



3 j _j_ v-j SjSj 



J/3 W;<4 



54 



m = 3 



(0|3) 



DV 



t=0 



oo 2k 2 . „ , 

+ E —kT^ 99 ^ [ Fli?2 " ^1(^)^1] + e E {{cdd) k Ff - 3F 1 (cdd) k F{ 



F 3 -(Sd)F 2 + -(Sd) 2 F 1 + 



k>l 



6 ^2 fc! 



E 



s 



2t 



6 
+ 



E 



5? / 2a? 6 



(Sd)£ 



2 [ ^2 1 l 3 

Mi Mf 



2Mi 



k k 



V 



Mi 



+ 









6 


r if)' 





3 2 



E 



2SfS J 1 a\ 



Mi \ 2a 2 



+ 

ij a ij 



2SiSjSk J 



2 \ -» 5f 

3 ^ M, 



Ei + 7 E 



at 4 ' aija ik 



j<k 



J 



+ E s i s j 



+4 V SiSjSk I — — !- 



1 1 

+ TT + 



MiOijUik MjCtjiOjk M k akiakj J 

(III.3.10) 

As conjectured by R.Dijkgraaf and C.Vafa [25], this expression coincides (up to rescalings S \tfi\ = 2S, 
= 2.F) with the result of calculation hi the framework of Seiberg-Witten theory |35|-|39|. 14 



4 According to |51|. 



T 



(0|3) 



i i^j i<.j<.k 

s / 



Mi = 



V ^ 1 1 X ^ 1 \ _ _\_ 



x - 1 5 

> -r + 7 



E — 



(III.3.11) 



+ 



E- 



111 
+ T7 h 



MiOLijOLik MjajiOjk MkCtkiCtkj 
and the two expressions for Mi coincide because of the sum rule |39|. 



El _ 1 1 1 1 

a%Mj ~ Mi ^ a?,- Mi ^ Ctii Ctik 



jjii *J J j^i 13 

obtained by summation over all the residues a* of the integral 



j<k 
j.k^i 



dx 



dz 



( 



E^ + E- 

— ' rvr. — < rtiAi 



OtijOtik 



The identity between 1111.3.1011 and 1111. 3. lit was first demonstrated by V.Pestun (January 2003, unpublished). 
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3.2 Example of n = 2 (cubic potential W(z), 3T 3 = «, M x = -M 2 = Ka i2 ) 



In the case of n = 2, the only non- vanishing cr-parameter is a\ 
formulas we easily reproduce the old results [2*3] : 



(i) 



-a 



(2) 



(0|1) 



DV 



t=0 



-5iW(ai) - S 2 W(a 2 ) 



(ill) 



t=0 



6kq 3 2 



r (2|l) 



(Si-S 2 )«3 6 » (2) _ 35(5i-5 2 ) 



t=o 



6! • 3 6 K 3 a? 2 



6k 3 a 9 



12 



F. 



(0|2) 



DV 



(1|2) 
DV 



t=0 



t=0 



--(5 1 2 -4S' 1 5 2 + S' 2 2 )loga 12 
(S? + S 2 2 )«3 4 » (1) 



4! • 3 4 • K 2 a 



12 



2S 1 S 2 /^ ((4 1 )) (1) + 2((3 2 )) (1) + ((2,3 2 )) (1) + 2 ((l 1 , 3 3 )) (1) 



2!-3 2 



3!-3 3 



2^,6_ 



12 



Z4536 
U-3 5 



(S 2 + S 2 2 ) 



1 



2™6_ 



re* a 



12 



9 

-(Sf + S|) — —51^ 



18 



3 4 



l/a±2, and from above 

(III.3.12) 
(III.3.13) 
(III.3.14) 
(III.3.15) 



(o) 



, 4-3 2-18 4-162 . . c , 

1 + — + —T7T- + ~A + 2 S X S 2 



(0|3) 



DV 



t=o 3 2 



6 Sf — Sf 5SiS 2 (Si — 5 2 ) 2 / 3 3 15 c c , , 

— ^ " " 3^4 I 51 " 52 " T 1 2 ^ 1 " 52 - 



12 



KCt 



12 



(III.3.16) 



(III.3.17) 



56 



Part IV 

Comments 



1 Generalized sum rule ( JIL1.25f ): derivation and applications 

The generalized sum rule (|ll.l,25j) can be derived in many different ways (see, for example, |41j for 
various derivations). Here we present yet another possible derivation. For the sake of simplicity, we 
put in this section the coefficient M equal to 1. The dependence on M can be easily restored. 
From the Toda chain equation ()L2.13j) . one can obtain the recursion relation for the density 

p N+1 (x) + p N ^(x) - 2p N {x) = (IV.1.1) 

All the densities with N > 1 can be found directly from this equation, but we use another approach. 
That is, we introduce the generation function for the Gaussian densities, 

oo 

h(x) = ]T PN (x)X N (IV.1.2) 

N=l 

It is easy to see that equation (|IV.1.1|) leads to the differential equation for this generation function. 
More precisely, one gets 

2sinh 2 (<^)^J e 2 V (x) =d 2 x J e 2 V {x) (IV.1.3) 

where 

J e2v (x) =sinh 2 (^)/ e 2„(x) (TV. 1.4) 
Making the inverse Laplace transformation for J, one gets 

2sinh 2 (^)^J(p) = p 2 J(p) (IV.1.5) 

where J(p) is the inverse Laplace image 

^ POO ^ 

J( x ) = / e~ xp J(p)dp (IV.1.6) 
Jo 

This differential equation can be solved, and one gets the general solution 

WP)= P . h2( \ } f( P ) (IV.1.7) 
sinn [ip) 

To find the integration constant f(p), one can determine the generation function for N = 1 and 
compare it with the explicitly calculation. This gives f(p) = \. To get the generalized sum rule 
(jll,1.25|) . one should consider the Laplace image of I\, 



(1 ~A)^ 
i.e. the integral 

e -P* e -— (IV.1.9) 



(1-A) 2 

which is divergent. To deal with this integral, one should treat it as a formal series in 1/x. Then, 
after simple calculations, one gets 



T , s 1 x\ f 1 

hix) = v^a^w J 

i ^ n^ + t^ ~_ £ * xN (IV - L10) 

—p: k —lie 2 — 

x 2 -t 2 i ' 2t 2 
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27T N=l 



X 
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and the coefficients of expansion in A give the desired sum rule. Unfortunately, we know any 
generalizations of the derivation neither for the Gaussian multidensities, nor for the non-Gaussian 
(multi) densities. 

Using the sum rule, one can easily restore the genus expansion. Introducing the variable x such 
that x 2 = jj, one gets the genus expansion as the l/7V 2 -expansion 



Pn{x) 



(I 



: 2_ t 2 



N 



1 e~ 2- 



dt 
2t? 



e^exp (JV^UaJ 
fe=i 

2t 2 oo £2 

exp(]T (~2 " X ) 



x t 2 \ 2fc+l 2 



(2fc + r 



2fc+l 



fc=l 



xV x 2 — 4 



2XV27T 



ex p(H 

fe=l 



t z \2fe+l 



52 _ 4 



J 



(2k + l)N 2k 
2 



1 e~T 



e 2 



2*2 

dt 
2t2 



(2fc + l)iV 2fc 



e x ~ 4 e 2 



_t 2 



(IV.1.11) 



Expanding this formula into the l/iV 2 -series gives the genus expansion, and one can get all the ex- 
pressions from table (|1.2|) . 

One can find the expression for the leading term of the polynomial Q in (jlll. 1.9f) . For genus k, one 

has 



NxVx 2 



(4fe-l)!! 



(4fc + l)!! A 
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1 (2i + l)(2Jfe-2* + l) 



iVx /(4fc - l)!ffi 2fc - 2 



6A- 
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+ (4A: + 1)!!^ 



J (2i + l)(2k -2i + 1) 



+ 0(x 



2Ai-6 



where y 2 = x 2 — 4. 
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Evaluation of multidensities in Gaussian model 



In principle, all the multidensities in the Gaussian model can be derived step-by-step from the Virasoro 
constraints, as are the first few in s ll.4l But it is desirable to get expressions like (|II.1.25j) which unify 
all genera. Here we consider only the two-point function. 

• There is a sum rule for the two-point function for N = 1 



p ( - 2 Xz 1 ,z 2 )4 1 z% 2 dz 1 dz 2 = 
(2fc 1 + 2fc 2 -l)!! k u k 2 odd, 

(2h + 2k 2 - 1)!! - {2h - 1)!! - (2k 2 - 1)!! ki, k 2 even, 
else 

The Virasoro constraints for the Gaussian potential can be combined as follows 

(xV x - d )z N = V 2 X Z N 



(IV.2.1) 



(IV.2.2) 
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Thus, for the disconnected two-point resolvent, one gets 



p {2) (x,x) + p (1) (x) 2 



x z + t z 



,2 _ t 2 



* 1 nnT t 2 \ _£ dt 
x 2 ' 2t 2 
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• From the Toda chain equation (jl.2.13f) . it is easy to derive the following equation for the two-point 
resolvent 

(d x + d y fp%\x,y) d 2 xP %\x)d 2 y p^\y) 



p$ +1 (x,y) + p)^ , __ 1 (x,y) - 2p^'(x,y) 



S 2 ) 



(2), 



(IV.2.4) 



N N 2 

This equation is consistent with known low genus densities, that is, they satisfy the equation in the 
leading order in 1/N. 

• Using the harmonic oscillator calculation as in |41j . one gets 



TTN ] {y, z ) = w(y + z) - G N (y,z) 



where 



N-l N-l 



Gn{v,z) = exp( )y B(n,n;y)B(n,n;z) + 2 ^ ^ B ( n i s 5 v) B ( n i s 5 

~ n=0 n=0 s=n+l 



(IV.2.5) 



(IV.2.6) 



Here B(n, s; x) is some (Hypergeometric) function, a proper matrix element for the harmonic oscillator 

r n+s— 2i 

(IV.2.7) 



B(n, s; x) = V nisi 



i=0 



i\(n — i)\(s — i)l 



Using the Laplace transformation, one can get the two-point resolvent 

roo pod 



P (2) (s,t) 



< ys+zt ^ (2 \y,z)dydz 



'0 JO 

• Using the integral representation of the Hermite polynomials H n 



H n (x) 



(-l) n n 
2m 



dz 

yU + l' 



-z 2 /2-zx 



one can represent the function Gtf(y,z) as the four-fold contour integral 



G N (y, z) = exp ( - ^ - ) ( I dzf ^ 
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(IV.2.c 



(IV.2. 10) 



For the generation function, one gets the expression 



if (A) = exp ( - V -±— ) J] Gn^ N = ( y ^) 4 ^ ^ dpidpae 
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The full generation function for the two-point correlation function is 



V n {2) X N - 



A (« + z) 2 l + A y 2 +z 2 

2 e 2 !- A - e 2 if (A) 
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(IV.2.12) 

It is easy to check that this formula gives an answer consistent with (|IV.2.6|) for N from 1 to 3. 

This formula is much more involved than the corresponding expression (|IV.1.8|) for the one-point 
densities, and this is not immediate to get the genus expansion for the two-point function from it. 
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3 Continuous limits 



There is an infinite variety of continuum limits, and only very few have been partly investigated. 
Among these are: 

• Naive large- ./V (planar or t'Hooft') limit (541 155j 

• Double-scaling limit jSHl EH EEl 

• BMN limit 

• DV limit Hi Hi S3 

• ... 

In all these limits the FSHMM converts into some other matrix model. Thus it is natural first to 
study these other models by themselves and only after that perform a detailed and justified analysis 
of the variety of continuum limits. This is a subject of separate publications. 

Also, particular representatives of these other matrix model families are identically equivalent to 
FSHMM, what allows to analyze the limiting procedures inside these families: various continuum 
limits are associated with different generalizations of the FSHMM. One of such examples is provided 
by the family of Generalized Kontsevich models |14| Ej], in which the Gaussian Kontsevich model 
coincides with the FSHMM. 
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